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The Voyager Plasma Science (PLS) instrument is capable of returning integral (DC) current measurements as 
well as normal plasma measurements. These measurements are similar in some respects to measurements 
made with a Langmuir probe or a retarding potential analyzer, although there are significant differences. The 
integral measurements have been made during calibration sequence in the solar wind, during Cruise Science 
Maneuvers, and within the magnetospheres of Jupiter and Saturn by Voyager 1. After the failure of the PLS 
experiment after the Saturn encounter, that instrument was placed in the DC return mode and returned 
possibly usable data from early 1981 through early 1985. The DC return measurements are difficult to 
interpret and are above threshold values only for relatively large fluxes; the determination of the measured 
current level is dependent on the operating temperature of the preamplifiers which further complicates the 
interpretation. Nevertheless, these measurements can be used to determine the efficiency of the suppressor 
grid at preventing the loss of secondary electrons off of the collector plate. Some DC return measurements 
have been invaluable in aiding in the interpretation of some electron plasma measurements which had not 
previously been understood. We find that the electron spectra can be significantly modified by the presence 
of second generation secondary electrons produced by the production of either first generation secondaries or 
photoelectrons on the support ring of the negative high voltage modulator grid within the instrument housing. 
This accounts for electron spectra obtained when the PLS side sensor is illuminated by sunlight while in one 
of the electron modes. The effect also explains the characteristics of electron spectra obtained by Voyager 2 
in the plasma voids near Ganymede during the Jupiter encounter and in the "charging region" while that 
spacecraft was outbound from Uranus. As a result, we have also identified other regions in which the 
Voyager spacecraft was charged to —1 kV and have noted regions in which care must be taken in the 
analysis of the supra thermal component of the measured electron distributions. 


Introduction 

The Plasma Science (PLS) experiment on the two Voyager 
spacecraft has returned a wealth of scientific information 
concerning the state of the interplanetary medium [Gazis, 1984] 
and the plasma populations of the magnetospheres of the outer 
planets visited to date [ Belcher , 1983 and references therein; 
Sittler el al., 1981; Lazarus and McNutt , 1983; Sutler et al. t 
1983; Richardson, 198 ; McNutt et al. t 1987; Selesnick and 
McNutt , 1987]. In addition to the plasma measurements made 
most of the time by the instrument upon which all previously 
reported analyses are based, there are integral (DC return) 
measurements. These measurements effectively bypass the 
action of the modulator grid in the instrument and are only 
taken over a limited dynamical range [Bridge et al., 1977]. 
These measurements are also difficult to interpret as is made 
clear in the body of this paper. However, they have proven to 
be an invaluable diagnostic tool in better understanding the 
effects due to the production of secondary electrons within the 
instrument housing itself. In particular, there are electron 
spectra which do not "make sense” given the normal operating 
characteristics and environment of the PLS instrument Using 
the DC return measurements and the production of 
photoelectrons within the instrument as guides and diagnostics, 
we have been able to explain semi-quantitatively data which had 
previously defied explanation. 


j Copyright 1988 by R. L. McNutt, Jr. 


DC Return Measurements 

DC return currents and normal plasma data are returned from 
the A and D cups of the PLS instrument in the DC return mode; 
plasma data from the B and C cups is replaced by A cup DC 
return and D cup DC return currents, respectively. The 
quantization steps for DC return currents are large to 
accommodate both large fluxes and currents of both positive and 
negative sign. The nominal sensitivity/range of the DC 
amplifiers (there are 2 DC measurement chains) is ± (1C 10 to 
10‘ 7 ) amperes [Bridge et al ., 1977]. In the main PLS analysis 
software package (VGRANL), there is a conversion lookup 
table from DN to femtoamperes. This table, in the common 
block PLCONS, has negative currents from - 1.00 x 10 5 fA 
(femtoamperes; 1 fA - 10' 15 A) for DN*0 to - 1.00 x 10 8 fA 
for DN»127 and positive currents from +7.5 x 10 4 fA for 
DN*128 to +1.00 x 10* fA for DN-255. A listing of these 
currents is included as Table 1 and a graph of the absolute value 
of these currents versus DN is included as Figure 1. The 
calibration curve from which this table was constructed is not 
known, but it is probably based upon a laboratory calibration 
curve (A. Mavretic, private communication) which was made at 
a temperature of +25° C (this is the most detailed curve which 
is available to me and is for SN003, the flight spare; there is 
also a '’quick check” measurement which was made for SN002, 
launched on Voyager; SN001 is on Voyager 2 - the spacecraft 
numbers and instrument serial numbers really are reversed). 
This curve is appended as Figure 2; note that for negative 
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currents there is a small difference between the conversions for 
the A and D cups. 

From another laboratory notebook (A. Mavretic, private 
communication), a full set of Ipc calibration have been found 
for SN002 (also designated Flight 1 or FI). These are from the 
summer of 1977. Figures 3 through 5 plot DC current (positive 
and negative) versus DN for both A and D cups. These were 
made for three different temperatures: -10°C, +27°C, and 
+50°D. Figures 6 through 9 show DC current versus 
temperature for constant DN values. Figures 6 and 7 are for the 
A cup (negative and positive current) and Figures 8 and 9 are 
for the D cup (negative and positive current). The calibration 
changes as a function of temperature (the preamplifier 
temperature) can be significant for the DC measurements; this is 
(thankfully!) not the case for the AC, i.e., plasma, 
measurements. 

DC return spectra typically show a DN value of 128 (positive 
threshold) in the D cup and a flat spectrum at the upper positive 
levels in the A cup (from DC return measurements in the L and 
M - positive modulator voltage - modes). At lower positive 
voltages in the A cup, higher positive currents indicative of the 
solar wind are apparent (at small enough heliographic 
distances). The currents at large potentials in the A cup have 
been ascribed to photoelectron production in the A cup, 
presumably off of the gold-plated collector plate. The D cup is 
usually in shadow so no photoelectrons should be produced. 
Considerable structure is found in the DC return spectra (all 
modes) from the Voyager 1 encounters with Jupiter and Saturn 
when particle fluxes into the cups were large; the PLS-A science 
link (the nomenclature may have been different for the Saturn 
encounter) included DC return measurements about every 45 
minutes. No DC measure menu were made routinely during the 
Voyager 2 encounters with Jupiter, Saturn or Uranus. The few 
DC return measurements made during the Voyager 2 encounter 
with Jupiter occurred far from the planet when particle fluxes 
were low and exhibit the same flat structure as found in the 
solar wind. 

R. S. Selesnick attempted to verify the photoelectron 
hypothesis by plotting the DC return current from the A cup 
(channel 128) as a function of spacecraft distance from the sun. 
At these voltages, there should be no contribution from the solar 
wind, even if the particle flux in the wind is very large. 
Selesnick's plou are appended as Figures 10 (Voyager 1) and 
11 (Voyager 2). Although the current tends to follow an r~ 
line for both spacecraft, there are unexplained deviations. Two 
possible effects not taken into account are i) the angle at which 
the sunlight is incident upon the collector and ii) changes in the 
DN to current conversion as the spacecraft moves out from the 
sun and cools (the calibration curve for a temperature of +25° C 
is used to do all of the conversions). The last point on the 
Voyager 1 plot is for day 236 of 1984 (corresponding to a 
heliocentric distance of 20.8 AU). In spot checking the printout 
of currents made at JPL (O. Divers, private communication), I 
found that in the L mode measurements in the A cup, the DN 
values in channel 16 are down to 128 by this time (Selesnick 
was using M mode spectra). In the lower channels, DC levels 
apparently due to the solar wind are being measured. The DN 
value in channel 1 is down to 129 (1 DN above positive 


threshold) on day 25 of 1985. The corresponding heliocentric 
distance is 22 AU. Information on the solar wind velocity 
and flux (hence, density) Is apparently available from the 
DC return measurements through this time. The M mode 
measurements seem to be clamped to a constant value indicative 
of the solar wind; in other words, the M mode measurements do 
seem to have been effected adversely by the failure. 

[n Voyager Memorandum #153, I gave a rough estimate that 
the DC return values should be near threshold at -24 AU. This 
distance was reached near the end of 1985, so the available 
numbers are roughly consistent At 400 km s’ 1 , 10 5 fA 
corresponds to a charge density of —0.23 cm ' (see Table 3 of 
Voyager Memorandum #99). 

Photoelectron Yields 

The photoelectron yield can vary as sec0 although this effect 
is depend*"* on both polarization and wavelength of the incident 
light in a complicated fashion (see Whipple [1981]). The 
energy distribution of photoelectrons can be well represented by 
a Maxwellian energy distribution, 

/(£) = £Xa e * e ° (l) 


Here J/(E) = 1 and £o is the "temperature" of the 
photoelectrons. 

In boating with both photoelectron and secondary electron 
emission, a flux "distribution function" is used [Grard, 1973; 
Hachenberg and Brauer, 1959]. In this case the Maxwellian 
energy distribution is 

j(E) = E/E& e' B£ ° ( 2 > 


Again, the normalization is such that 



1. Note that the 


most probable energy - defined at the maximum of j{E) ~ is 
E 0 . Hence, when a "most probable" secondary electron energy 
is quoted in the literature, it is almost always with respect to 
;(E) and. therefore, is the same as the secondary "temperature." 

Grard [1973] finds values of £ 0 ("most probable energy") of 
0.88 eV for gold (appropriate for the collector plate) and 1.02 
for graphite (appropriate more or less for the black conducting 
paint - up to 12% carbon by weight (A. C. Whittlesey, private 
communication, Feb. 20, 1987) - used on the Voyager 
spacecraft; secondary and photoemission characteristics of the 


paint are not known). 


PL S Instrument Temperatures 

Three instrument temperatures are routinely telemetered to 
Earth in the PLS engineering records. Knowledge of the 
temperature decreases during the missions of Voyagers 1 and 2 
as they move away from the sun is of both general interest as 
well as necessary for the proper interpretation of DC return 
measurements: the calibration curves for translating DN (data 
numbers) to Idc values is dependent upon the temperature of the 
preamplifier and measurement chain. 

The temperature readouts from the PLS instruments on both 
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spacecraft have been indicating decreasing temperature with 
increasing distance from the sun. Temperatures and other 
readouts are returned in engineering records and are plotted to 
produce engineering plots. Such plots have been made for both 
spacecraft for a great deal of the two missions and are in a 
notebook in the programmers’ office. Printout of the 
engineering records consists of 9 integer fields (IENG) and 3 
real fields (TEMPS). The integer fields are data numbers, and 8 
of these are plotted on engineering plots with the not terribly 
illuminating labels ESUP, VREF, VCOM, PSUP, TSEN, 
TPMC, TMOD, and MODE. The non-temperature data is part 
of the PLS analog multiplexer output. The latter is comprised 
of the PLS backup science data ( M a 0 to 3 volt representation of 
the four decade plasma measurement to be used in the event 
that the PLS ADC ceases to function'’), the internal reference 
voltage of 10 volts, the suppressor voltages, and a combination 
of the positive and negative 15 volt supply voltages; 
correspondence with the engineering plots is obvious given this 
information. Backup science data DNs are not routinely plotted; 
backup science data is also not well documented at this 
juncture. What documentation I have found is in JPL document 
618-505, Vol. IV, Rev. B. 

The temperature data are interpreted as follows 


Label 

Description 

IENG( ) 

TEMP( ) 

TMOD 

Modulator Temperature 

1 

1 

TPMC 

Preamp and Measurement Chain 

2 

2 

TSEN 

Sensor Temp. (A cup Housing 7 ) 

3 

3 


The units of TEMP are °C; the conversion algorithm from DN 
to °C is not documented but buried in the subroutine MJSTEM. 
The DN values are first converted to resistance values. These 
are then converted to temperatures via interpolation from a 
lookup table stored with the subroutine. The sensor temperature 
conversion is different from the other two. Lookup tables and 
graphs for the conversion for the sensors on S/C 31 are included 
as Appendix A (from JPL document 618-301), 

Engineering plots were not routinely made during the 
encounter periods; they have been recently run for one day 
periods during the encounters. Variation during these periods 
was minimal; typical temperatures from the encounters are listed 
in the following table. 


S/C 

Due 

Planet 

TSEN TO 

TPMC CO 

TMOD (*Q 

i 

1979/64 

Jupiter 

•123-58 

9.07 

11.18 

i 

1980/318 

Saturn 

• 148-22 

4.14 

6.95 

2 

1979/189 

Jupiter 

•12X39 

9.07 

11.18 

2 

1981/237 

Saturn 

-147.17 

4.84 

6.95 

2 

1986/ 24 

Uramu 

-16035 

3.43 

65)5 


Selesnick's numbers (referred to above) have not been corrected 
for the temperature shift; how much this will "help" his curves 
is not yet known. 


Photocurrent Loss From DC Return Measurements 

At the Earth Garrett [1981] lists the photoelectron flux from 
solar spectrum photons incident on gold as 
/p* = 2.9 x 10- 9 A cm" 2 (see his Table 4). The effective area of 
each of the main sensor cups is A^ = 66 cm 2 (see p.12 of 
Barnett, [1984]), and the photoelectron flux at Jupiter is down 
from that at Earth by a factor of 5.2 -2 . This should produce a 
photocunent at the collector plate of 

Ipk - jpk&wtk r~ 2 (3) 

or 7.08 x 10 6 fA. I have looked at the L mode DC return 
measurements from the Voyager 1 Jupiter encounter (DC return 
measurements were made every 45 minutes as part of the PLS- 
A science link). The following table lists DC return times and 
the A cup current and the D cup current in channel 16 (the D 
cup current was at the positive threshold of 7.5 x 10 4 fA for 
most of the spectra). Note that the currents listed are nominal 
values using the lookup table exhibited in Table 1 and Figure 1. 


Time 

U.16 fA 

DN, A 

^D,16 fA 

62 0601 

1.71 x 10 s 

192 

7.50 x 10 4 

63 1924 

1.38 x 10 s 

191 

7.50 x 10 4 

64 0503 

820 x 10 s 

184 

7.50 x 10 4 

64 0737 

1.38 x 10 s 

191 

7.50 x 10 4 

64 0921 

1.38 x 10 6 

191 

7.50 x 10 4 

64 1426 

1.38 x 10 s 

191 

7.50 x 10 4 

64 1829 

-2.40 x 10 s 

66 

-2.0 x 10* 

64 1922 

-2.40 x 10* 

66 

-4.9 x 10* 

65 0629 

4.14 x 10 5 

164 

-2.0 x 10* 


The difference between the fluxes into A and D should be 
representative of the current measured due to the loss of 
photoelectrons within the instrument (i.e. those not returned to 
the collector plate by the potential on the suppressor grid). If 
we take 1.5 x 10* as a representative photocurrent, this implies 
an internal loss of ~ 20 %, consistent with more detailed 
calculations on secondary electron losses in the side sensor. 
From Voyager Bulletin #37, entrance into solar occultation was 
1637 and exit from solar occultation was 1849. This is 
consistent with the no photoelectron production and the 
measured negative current at 1829. The very large negative 
currents at 1922 of day 64 are apparently due to photoelectron 
production in the side sensor; due to a maneuver, the angle 
between the D cup normal and the sun decreased to just under a 
critical angle discussed elsewhere. The large negative currents 
on 65 are presumably from the hot electron population near the 
spacecraft, but this needs further investigation (and use of the 
proper DN to current curve). Four DC return measurements at 
Saturn on day 317 of 1980 (at 0736, 0738, 0828, and 0919) all 
gave the D cup at positive threshold and the A cup at 
5.80 x 10 5 fA, somewhat larger than predicted than a simple r' 2 
scaling from Jupiter would predict, but, again, this could be 
related to the change in the preamplifier temperature. 
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There are a few other things to note about the Voyager 1 
encounter with Saturn. The first is that there was one 360° yaw 
turn at Saturn, similar to the PLOMAN and FLOMAP 
maneuvers by Voyager 1 at Jupiter. For the record, appropriate 
information about the maneuver (FSMAN-5) is as follows. 


Voyager 1 Yaw Maneuver at Saturn 

Type of Rotation 

Start Time of Rotation 

Roll 

318/2117:02 

Yaw 

318/2128:24 

Yaw 

318/2146:29 

Roll 

318/2210:42 


My notes indicate that the first yaw went through 132° and the 
second completed the 360° turn. 

In addition, it should be noted that the sun occultation ingress 
for Voyager 1 was 1980-318/0153:23 and the sun occultation 
egress was 1980-318/0240:01. In occultation, there is no 
evidence from a cursory examination of the 1800 plots that any 
charging of the spacecraft to high negative levels was occurring. 
This is probably not surprising in that the bulk of the electrons 
had a relatively low (< 50 eV) temperature during this time. It 
should also be noted that of DC return spectra were acquired 
around - 318/0200 within the occultation period. The spectra 
from the A cup do show that the photoelectron current from the 
collector plate was turned off. They also indicate large negative 
currents at positive potentials up to several tens of volts. If 
such a signature has occurred elsewhere in the DC return 
spectra, I have not yet run across it 

Secondary Electron Emission From the Collector Plate 

I am also investigating the possibility of "experimentally" 
obtaining the suppressor efficiency from looking at DC return 
measurements made early in the mission when sunlight 
ill uminatMi the D cup collector plate. In addition, George and I 
are working on a quantitative model of secondary electron loss 
from the D cup. Secondary electrons can also be modeled with 
Maxwellian distribution functions [Whipple, 1981]. 

The collector plate is electroplated gold on an electrolous Ni- 
Cu phosphate substrate on Al. The substrate prevents the 
"leeching" of the magnesium structure into the gold plating. 
Plating details (from the instrument specifications) are as 
follows: 1) flash Cu to thickness of 100 ±50 micro in., 2) 
underplate non- magnetic electroless Ni per Mil-C-26074 Class I 
to thic k"*” of 500 ± 100 micro in., 3) overplate bright Au per 
Mil-G-45204 Type I Class II to thickness of 100 ± 25 micro in. 
(from drawing 21-8001, rev. C). Hence the gold plating has a 
thickness of at least 2 x 10~* cm. 

Sec o ndary electrons are emitted from a surface layer much 
thintwr than this, - lO^cm (see p. 423 of Dekker [1957]); hence 
the secondary emission characteristics of gold are appropriate 
for determining the secondary yield and energy spectrum. 


The most extensive compilation of the secondary emission 
characteristics of iw etaia appears to be due to Kollcuh [1947, 
1956]. The best "temperature" for secondaries emitted from 
gold is -2 to 3 eV (from table 1 of [ KollaiK 1956]) obtained 
from a primary electron energy of 150 eV. The secondary 
emission spectra are not very sensitive to the energy of the 
p rimar y for primary energies above - 20 eV [Hachenberg and 
Brauer ; 1959] and at least up to - 1.6 MeV [Schultz and 
Pomerantz , 1963]; these are well outside the energy range for 
which the total secondary electron yield is significant (see 
below). A temperature of - 3 eV is also consistent with the 
theoretical treatment of Chung and Everhart [1973]. For 
quantitative modeling of the loss of secondaries from the 
collector, we have assumed a temperature of 2.5 eV; as this is 
small compared to the instrument suppressor voltages of - -8 V 
and -95 V, small errors in the temperature have a negligible 
effect on our results. 

The loss of secondary electrons (i.e. those produced at the 
collector and not returned to the collector by the suppressor 
grid) was considered in detail for die PLS instrument by Sittler 
in 1979 (E. C. Sittler, Jr., private communication). The 
geometrical problem is illustrated in Figure 12 (figure due to E. 
C. Sittler). For a given radial location r on the collector plate 
(of radius r c ), an impact by a primary electron will cause the 
emission of secondary electrons with an energy distribution as 
discussed above and an angular distribution which is isotropic 
[Whipple and Parker , 1969] (for an extended emitting surface) 
with respect incidence angle 0 and azimuthal angle +] (early 
work, see e.g. Hachenberg and Brauer [1959] and McDaniel 
[1964], suggested a cos 0 emission distribution, but isotropy 
seems to now be accepted by most workers in the field [ 
]). Consider a secondary electron of energy E t emitted at angle 
0. The trajectory will follow a straight line in the field-free 
region between the collector plate and the grounded grid (Grid 1 
- this is the "normal suppressor" electrical configuration). Grid 
2, held at the suppressor potential V^p will decelerate the 
electron motion perpendicular to the collector plate, resulting in 
the parabolic trajectory shown in the figure for Ep € < 
higher energy elections penetrate the potential barrier and are 
i mme diately lost from the system. The parabolic trajectory will 
again intersect grid 1 a distance R from the position at which it 
first crossed the grid plane. If d ltv is the distance between the 
grounded grid and the suppressor grid (equivalent to z 2 - z i in 
Fig. 12), then an elementary calculation shows that 

R = 2 d^, ■ sin 20 (4) 

The electron will again move in a straight line until it crosses 
the plane of the collector plate. If d (z x in the figure) is the 
separation of the collector plate and the ground grid, the 
electron will intersect the plane of the collector plate a distance 
a (equal to Ax in the figure) from where it was created where 

a = 2d tan 0 + R (5) 

If the impact point is outside of the collector, the secondary 
electron is lost, giving rise to a spurious (positive) current in the 
data. From Figure 12, it is obvious that the electrons emitted in 
a cone of half angle 0 map to a circle in the plane of the 
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collector plate. As long as a- r <r e there will be a range of 
azimuthal angles <j> such that |0| £ 0n»x for which the 
secondaries will be lost For larger values of a (corresponding 
to smaller values of 0 ), all of the secondaries will be lost 
Noting again that both the angular and energy distributions of 
the secondary electrons is independent of the energy and angular 
distributions of the primary electrons, we can calculate a solid 
angle corresponding to secondaries which are lost 

Define the angular limit G^E,) as the smallest value of 0 for 
which a = r c using equations 4 and 5 (found numerically); 
define 0 2 (E,) as the largest value of 0 for which a = r c or it/2, 
whichever is smaller (N.B. 0 considered as a function of R is a 
double-valued function). Define 4>o by 


MG.E,) 3 cos 1 


r 2 - r 2 — a 2 
2 7a 


for a < r + r e (6) 


(0.E,) m n for a £ r + r c (7) 

The angular l imi ts are, of course, dependent upon the energy of 
the secondaries due to the parabolic section of their trajectories. 
The solid angle for lost current (integrated over the secondary 
electron energy spectrum) is 

-V W 

0^) = 2 f AE,)dE, f m MO.EJ sin0 d0 (8) 
i e/irp 

where A(0) is the angular distribution of secondaries, and j(E j) is 
their energy distribution function with respect to current, given 
by equation 2. We have computed ft^r) both for the case of 
isotropically emitted secondaries, i.e., A(0) = 1 and for a cosine 
distribution (normalized to unity over a hemisphere), 
/t(0) = 2cos0 for the various operating modes of the instrument 
There are four relevant operating modes of the instrument for 
which fWr) has been evaluated: the normal suppressor 
configuration, described above, and the reversed suppressor' 
configuration, in which the grounded grid and suppressor grid 
are reversed; for each configuration there are two suppressor 
voltages, one for the low energy electron measurements (El) 
and one for the high energy election (E2) and ion (L, M) 
measurements (see Bridge et al . [1977]). Relevant parameters 
are listed in the table below. 


Mode 

Configuration 

V,* (V) 

d (cm) 

(cm) 

El 

Normal 

-8.05 

0.996 

0.417 

El 

Reversed 

-7.70 

0.000 

0.996 

E2.L.M 

Normal 

-95.2 

0.996 

0.417 

E2, L, M 

Reversed 

-85.9 

0.000 

0.996 


The distances between the grids are from Figure 2.12 of Barnett 
[1984]. The cups of the main sensor only have a "normal" 
configuration and use a suppressor voltage of -95.2 V (these 
numbers are from Voyager Memorandum #69 of 1/8/80 which 
indicates the data source as Mavretic’s notebook #4. The 
voltages quoted are within 50 mV for SN002 and SN003; it is 


assumed that SN001 has the same values, but these are not 
documented). 

Note that for $o = it, 0j = 0, and 0 2 = it/2, = 2n, i.e. 

the solid angle in a hemisphere. G. S. Gordon, Jr. has 
normalized fl^r) to 2n to find the fraction of secondaries lost 
per unit incident current as a function of radial distance r from 
the center of the D cup collector plate. E. C. Sittler, Jr. 
performed the same calculations in 1979; his results for the 
fraction of secondary electrons lost in the El and E2 modes are 
shown in Figures 13 and 14, respectively. Sittler normalized the 
losses to their maximum values and expressed radial distance as 
a fraction of the collector radius. The maximum loss occurs at 
the edge of the collector plate. In El Ed found a loss of 26.3% 
for the normal suppressor configuration and 12.5% for reversed 
suppressor configuration. For E2 the loss fractions at the edge 
are 25.3% and 3.5%, respectively. 

The results of Gordon's calculation are shown in Figures 15 
through 19. Again the normalized value of fijr) is plotted. 

In Figure 15, results for all of the cases noted in the table above 
are plotted. The notations "iso" and "cos" refer to isotropic and 
cosine angular distributions of the emitted secondaries, 
respectively. The notations "std" and "rev" refer to whether the 
suppressor voltage is on the standard or reversed grid. The 
suppressor potentials are also indicated. Figures 16 through 19 
show the loss fraction for cosine and isotropic secondary 
distributions for each of the four suppressor modes; the dotted 
lines indicate the isotropic results and the solid line cosine 
results. Note that in the standard configuration - denoted by "s" 
in the label on the graph - the isotropic distribution produces 
more losses as there are more electrons emitted at large angles 
of incidence which are not reflected back to the collector until 
they have moved relatively far from their point of emission. In 
the reversed configuration - denoted by an "r" - the cosine 
distribution gives rise to more losses as there are more "forward 
emitted" secondaries which simply penetrate the potential barrier 
in the forward direction; electrons emitted at large angles of 
incidence are immediately "shoved" back into the collector plate 
by the suppressor potential. Note that *T refers to the case of 
L, M, or E2 and "e" in the label refers to the El potential, e.g., 
the label "is" identifies the normal, or standard, suppressor 
configuration with -95.2 volts on the suppressor grid. 

Listings of the loss fraction as a function of radial distance 
across the collector plate are included in Appendix B. The file 
outrlmel5ci" contains numerical output for the isotropic cases; 
the file "outrlmel5cc" contains the the corresponding numbers 
for the case of a cosine angular distribution. The parameters 
used for each table are listed at the top: "ds", "d", and "volts’ 
correspond to d^ d, and respectively. The parameter 

M e0" refers to the assumed energy of the secondary electrons, 
here 2.5 eV for each case. 

We have not included either elastically scattered (reflected) 
primaries or inelastically scattered (rediffused or backscattered 
primaries) electrons [Hachenberg and Brauer , 1959; Dekker , 
1957]. Reflected primaries have the same energy as the incident 
electrons and display a maximum in the direction of the incident 
beam [Hachenberg and Brauer , 1959]. According to Whipple 
[1981], reflection is only significant at primary energies less 
than -10 eV. Backscattered electrons occupy the energy range 
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be tween -50 eV and the energy of the primary. They obey the 
same angular distribution law as the true secondaries 
[Hachenberg and Brauer ; 1959] and can contribute a net 
backscattered current - 0.3 of the incident current [ Whipple , 
1981]. As the true secondaries contribute most of the current at 
the collector, we only consider their effect here; it should be 
noted that a sizable fraction of the backscattered electrons will 
be lost in any event if the energy of the primaries substantially 
exceeds eV Si ^. 


Secondary Yield and Lost Secondary Current 


The loss of secondary electrons produced at the collector 
plate is measured as a positive current by the DC measurement 
chain, and, if modulated, as a positive current by the AC 
measurement chain (see Bridge et al . [1977]). This current is 
superimposed on the "direct” current produced by the primary 
flux of charged particles on the collector plate. It is convenient 
to consider separately the "direct'* or primary current as opposed 
to the secondary current This is done in terms of the secondary 
yield function 5 where 

U = (9) 


and the subscripts refer to the current density of the primary 
particles and that of the secondary particles [Hachenberg and 
Brauer, 1959]. In what follows we concentrate on the 
production of secondary electrons by primary electrons; at 
primary energies above several keV the production of secondary 
electrons by primary ions can also give rise to significant effects 
[Whipple, 1981]. 

The yield of secondary elections form primary electrons is, in 
general, a function of both the primary energy E p and angle of 
incidence 0,. However, the yield curve tends to be universal, 
i.e., independent of the emitting material, if parameterized in 
terms of the maximum yield and the energy at this maximum 
E m (see e.g., Dekker [1957]). The yield function is further 
characterized by the energies E f and E u corresponding, 
respectively, to the lowest and highest primary energy for which 
the secondary yield exceeds unity. For gold 6 m is 1.45 at E m of 
800 eV, and E t and E u are 150 eV and > 2000 eV, respectively 
[Kollath, 1956], An analytic form for the yield is given by 
Whipple [1981] which we use here 


5 ( Ep , B p ) = 
il - exp [- 


1.1148 m 

E m 

COS0, 

E 

» 

-2.28cos0 

* 

J l 


Em 


1 0.35 


11.3511 


(10) 


Note that for small energies or angles of incidence near 90°, 

5(£^6 P ) oc Ep. 

The response of the PLS side sensor to electrons is a 
convolution of the the distribution function of the electrons with 
the "response" function of the detector. The response function 
can be conveniently separated into a transparency function T 
and an area overlap function A; the latter contains the effects of 
the refraction of particle trajectories by the various grids as well 
as the decrease in effective collector area due to particle 


trajectories within the detector [Barnett, 1984; Barnett and 
Olbert, 1986]. 

If we define a coordinate system with the z axis along the 
symmetry axis of the side sensor, or D cup, then the angle of 
incidence of a primary election with respect to the collector 
plate is the same as that with respect to the entrance aperture of 
the detector (both the grid at the entrance aperture and the 
collector plate are at the spacecraft "ground" (which is the 
potential of the conducting surface of Voyager, the latter can, of 
course, have some floating potential with respect to the 
surrounding plasma). Hence 0, = tan’V/v,) where v, is the 
component of electron velocity into the sensor and v, is the 
component of velocity transverse to the sensor. The energy 
dependence of the yield follows from E p = in m* (v? + vj). 

For a given angle of incidence and given primary energy, the 
region of the collector plate illuminated by electrons passing 
through the aperture grid, will, in general, be the common 
overlap area of three circles whose centers are coline ar [Barnett, 
1984]. This is a two parameter function A{d mg d^ which can be 
described by the distance between the centers of the aperture 
image and the guard ring image d* and between the centers of 
the guard ring image and collector image respectively; the 
function can be written in closed analytic form [ Barnett , 1984]. 
The distances are offsets determined by the refraction effects 
which are, in turn, functions of v„ v„ and the 
modulator potential. The current measured due to electrons 
more energetic than eV mod is 

m ma mm 

Idc/k = e j dv.v, j j </*v,/(v) T A (11) 


where v m ■ (2eV^ /m t ) ia . The DC return current actually 
measured by the instrument is loc = in (//x>, + Ioc*» i), i.e. the 
average current in a given channel (when the modulator is on, 
which is usually the case). 

To calculate the effect of secondary current loss, we have 
modified the area overlap function as follows. We have 
numerically integrated the solid angle for secondary loss Q^r) 
over the overlap area to find an effective "loss area" A^ 

Aiau - J <*A (12) 

os*riapar + 1 

The secondary current lost due to primary electrons more 
energetic than eV^ is given by the current density at the 
collector plate multiplied by the secondary yield and the loss 
area, viz. 

/«« = « j dv, v, j j <Pv,f(y) T 5 Aun (13) 

— — 

(Note that the lost secondary current density per velocity space 
volume is simply the primary current density at the collector 
plate x 5 x 

In the DC return mode, the current measured for modulator 
voltage Vmai is given by loc - /**• In the normal plasma mode, 
the measured current is the modulated part of this current 

To incorporate this effect into the simulation of electron 


spectra, Gordon computed a set of eight loss tables (low *nd 
high suppressor potential; normal and reverse suppressor 
configuration; isotropic and cosine law for the angular 
distribution of the secondary electrons). The numerical 
integrations of eqn. 12 for these 8 cases were run on a Sun 
Microsystems Workstation on the weekend preceding May 18, 
1987. The jobs were run in the background beginning on the 
preceding Friday and completed late on Sunday, May 17. 

The loss tables have been incorporated into the PMODEL 
package via modified versions of the subroutine REDART 
(which reads the tables into the SUN from the disk where they 
are stored) and the functions ARS1DE and CPRESP. In 
CPRESP, the overlap area A (refer to eqn. 11 above) has been 
replaced by A - 8 A**, (refer to eqn. 12 above). The output 
from PMODEL then produces synthesized electron spectra 
which are corrected for the inefficiency of the suppressor. 
Listings of the relevant code are included as Appendix C. Note 
that the hand lettering on the listing for CPRESP is only 
applicable to using the code to simulate the effects of 
photoelectrons ; the code as printed on the laser printer was 
implemented for the case of secondary electrons made from 
primary electron (rather than photon) impacts. Equation 10 with 
b m - 1.45 and E m - 800 eV is used to compute the yield. 

Figure 20 displays a contour map of the normal area overlap 
function; the two axes are the relative displacement of the 
aperture and shield circle centers (DAS) and the relative 
displacement of the shield and collector circle centers (DSC) 
(these correspond to the coordinates X and Y, respectively, 
which appear in the call to ARSIDE in the function CPRESP). 
Figures 21 and 22 show the area (in cm 2 ) as a function of one 
of these two parameters (labeled M l" and "2") for a constant 
value of the other. The distance from one side of the aperture 
to the other side of the shield is 10.8 cm, and the distance from 
one side of the shield to the other side of the collector plate is 
11.5 cm; hence, these are the values of DAS and DSC for 
which the area overlap goes to zero. Parameter 1 is therefore 
DAS, and parameter 2 is DSC. A three-dimensional projection 
of the overlap area function is shown in Figure 23. 

Figures 24 through 39 are plots of A **, „ versus DAS (*'T ) 
and DSC (-“2"), in an analogous fashion to the plots of A in 
Figures 21 and 22. Each plot is labeled with a four character 
alphanumeric label: “vcan M , where M v" is "e M (El suppressor 
voltage) or “i" (L, M, and E2 suppressor voltage), "c” is V 
(normal suppressor configuration) or "r M (reversed suppressor 
configuration), "a" is *T (isotropic angular distribution of 
secondaries) or "c" (cosine angular distribution of secondaries) 
and "n” is "1" or “2“ depending upon whether DAS or DSC is 
the independent variable plotted along the abscissa. Note that in 
each case in which the independent parameter is DSC, the 
effective loss area is a minimum for DSC=0. As this parameter 
is increased, the losses increase (refer to Figure 15), but the 
overlap area decreases (Figures 21 and 22). Hence, there is a 
maximum loss for some nonzero value of DSC. Also note that 
the scale along the ordinate is different for Figures 24 through 
39. 

Figures 40 and 41 are contour plots for the E2, normal 
configuration case, for isotropic and cosine distributions of the 
secondaries, respectively. Corresponding three-dimensional 


representations for these two cases are shown (on the same 
scale) in Figures 42 and 43, respectively. Figure 44 reproduces 
Figure 42 but on the same scale as Figure 23; from this last 
comparison, it is obvious that the loss of secondary electrons 
from the collector plate can decrease the measured electron 
current but cannot change its sign (a seemingly trivial point 
which is important in accounting for some of the PLS electron 
spectra). 

We simulated electron spectra (E2 mode) by using PMODEL 
and including the loss portion for assumed isotropic and cosine 
distributions in the normal suppressor configuration. Figures 45, 
46, and 47 show simulated currents for an isotropic Maxwellian 
distribution of electrons at a temperature of 400 eV assuming, 
respectively, no losses, isotropic emission of secondaries, and a 
cosine law for the emission of secondaries. The losses decrease 
the measured currents by -20% or less. Note that the apparent 
discontinuity between channels 3 and 4 is due to the suppressor 
grid; the threshold of channel 4 is -92.3 volts whereas the 
suppressor potential is -95.2 volts. The currents in the first 
three channels are due entirely to modulation of the trajectories 
of the electrons in the input distribution with energies greater 
than 95 eV (this is the so-called feedthrough effect). To further 
investigate this effect, we also simulated a k distribution with a 
most probable speed corresponding to an energy of 9.85 eV. 
The case for no losses is shown in Figure 48. This spectrum is 
qualitatively similar to those observed in the vicinity of the 
“voids” in the Jovian magnetosphere by Voyager 2 [McNutt et 
al. 9 1987]. The high energy cutoff is due to numerical cutoffs 
used in simulating the spectrum. The spectrum including losses 
from isotropic secondaries is shown in Figure 49 and the loss 
current itself in Figure 50. 

Experimental Verification of the Suppressor Efficiency 
General Considerations 

The calculation of the apparent current due to the inefficiency 
of the suppressor grid is complex enough that it is desirable to 
verify it experimentally. Efforts are currently underway at 
M.I.T. to do so in a vacuum chamber using the side sensor of 
the Flight Spare PLS experiment (as of this writing, Sept 17, 
1988, these have not yet begun). In addition to what can be 
done in a laboratory (which can approach but not really 
duplicate the environment of interplanetary space), we can look 
for this effect in some of the plasma measurements themselves, 
in which the local environment is correct but the types of 
measurements were not really optimized for looking for this 
effect We noted previously that photoelectrons and secondary 
electrons have similar energy distributions for production from 
gold surfaces, the mean energy of the photoelectrons being 
smaller (0.88 eV versus ~2.5eV). These energies are small with 
respect to typical potentials within the PLS detectors; hence, the 
emitted electrons should follow similar trajectories. No 
experimental information on the angular distribution of emitted 
photoelectrons appears to be readily available; however, it is 
reasonable to assume that the angular distributions are also 
similar. We can thus use sunlight into the instrument hitting the 
collector plate as a probe to determine the suppressor loss as a 
function of the illuminated part of the collector plate. 
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When the instrument is in an election measurement mode, 
photoelectrons are emitted by the support ring for the negative 
high voltage grid whenever the angle between the cup normal 
and the sun is less than 78.5°. If sunlight also hits the collector 
plate, photoelectrons will also be emitted there, but the loss 
current is small compared to that produced by those 
photoelectron coming from the negative modulator support To 
eliminate this effect, the instrument must be in a positive ion 
detection mode. In this case, the only surfaces with negative 
potentials are the suppressor grid, which has a small cross 
sectional area, and the buckout grid. The latter has only a small 
negative potential which alternates with a positive potential. As 
the suppressor grid is between it and the collector plate, any 
electrons from the buckout grid cannot reach the collector plate 
in any case. 

The loss of photoelectron current from the collector plate 
cannot be affected by the modulator potential and will give rise 
to only a DC signal as the modulator steps. This is the source 
of the DC current in the high channels of the A cup measured 
during cruise. The geometry of the A cup is complex enough 
that a good mathematical model for the suppressor losses, while 
possible to construct, will also be very complex. Fortunately, 
there were a series of cruise maneuvers conducted during the 
cruise phase of both spacecraft during which the D cup collector 
plate is illuminated by the sun, and the instrument is in DC 
return mode. 

Cruise Maneuvers 

Several types of changes in spacecraft attitude are made 
during cruise. Early in the mission Cruise Science Maneuvers 
(CRSMVR) were performed, consisting of 25 roll turns 
(counterclockwise about the +Z axis of the spacecraft, which 
points opposite to the high gain antenna) and 10 yaw turns 
(counterclockwise about the +Y axis of the spacecraft, which is 
aligned with the projection of the science boom perpendicular to 
the +Z axis - see Figure 51). Later in the mission these were 
changed to mini-Cruise Science Maneuvers (Mini-CRSMVR) 
which included fewer revolutions of the spacecraft The 
original idea was to partially roll the spacecraft so as to bring 
the D cup directly across the sun during the yaws. As far as I 
can tell this was not done, probably as the result of operational 
constraints of which I am not aware. A listing of all maneuvers 
through a complete revolution for both spacecraft as well as a 
memo describing our original thoughts on PLS measurements 
(from J.D. Sullivan, dated June 14,1978) are included as 
Appendix D. 

I have looked in some detail at the CRSMVR yaw turns on 
Voyager 1 on day 258 of 1978. There is one sequence of DC 
return measurements (normal suppressor) for one full yaw early 
on day 258. Other measurements during the yaws are normal 
plasma and current calibrations (CURCAL) with both modulator 
on and off, all with normal suppressor. The minimum angle 
between the D cup normal and the sun is 43.3°, corresponding 
to the angle between the D cup normal and the +X axis of the 
spacecraft (refer to Figure 51). Sunlight reaches the collector 
plate in the D cup when the sun-cup angle is less than 63.1°. 
Due to the fairly rapid rotation rate of the spacecraft (-32 
minutes to rotate through 360°), this limits the available DC 


return spectra for which the theory and observations can be 
compared. 

The DC return measurements were removed from the 
sequence when the change was made to mini-CRSMVRs. From 
the maneuver list in Appendix D, this limits the comparison to 
11 maneuvers on Voyager 1 and 6 maneuvers on Voyager 2. 
The comparison is further limited by the distance between the 
spacecraft and the sun. The suppressor in the D cup is 
apparently more efficient at returning photoelectrons to the 
collector plate than is the suppressor in the A cup. In addition, 
the geometry during the yaws limits the illuminated area of the 
collector in the D cup. These effects combined with the 
relatively high threshold for detectable current in the DC return 
mode further limits the availability of data. As the spacecraft 
both receded from the sun, the observed current decreased to 
threshold before the CRSMVRs were changed to mini- 
CRSMVRs and the DC return measurements were deleted from 
the measurement sequence. The following tables list the 
CRSMVRs by date, the maximum DC return DN in the D cup 
in channel 16 of the L mode, and whether evidence of 
photoelectrons in the DC returns in the electron modes is 
present 


Voyager 1 

SCET 

Max DN (IDQ 

Photoelections 

Range (AU) 

1978-048 

172 


Yes 

2.4 

1978-257/8 

151 


Yes 

4.1 

1979-034 

146 


Yes 

5.1 

1979-288/9 

130 


Yes 

6.4 

1980-031 

N/A 


N/A 

7.1 

1980-168 

128 


Yes 

8.2 

1981-041 

128 


77 

9.9 


More information concerning the DC return levels after the 
Voyager 1 failure can be found in Voyager Memorandum #157. 
No maneuver data could be found from the tape logs on day 31 
of 1980. 


Voyager 2 

SCET 

Max DN (IDQ Photoelectrons 

Range (AU) 

1977-362 

177 

Yes 

1.9 

1978-275 

N/A 

N/A 

4.0 

1979-155 

136 

Yes 

5.2 

1980-024 

128 

Yes 

6.1 

1980-136 

128 

Yes 

6.7 

1981-148 

128 

Yes 

9.0 

1982-103 

128 

No 

10.0 


Data from a maneuver on day 275 of 1978 could not be located 
and the mini-CRSMVR on day 155 of 1979 does have DC 
return data during a yaw turn. The SOEs have not been 
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investigated closely to clear up these discrepancies. During the 
yaw turn on day 24 of 1980, the lower channels of the D cup 
registered 134 DN and fell off to 128 at the location of the solar 
wind peak in the normal plasma data for the D cup. The same 
behavior occurred on day 136 of 1980, with the first channel of 
the D cup showing a maximum of 132 DN. On day 148 of 
1981, the solar wind can no longer be seen in the DC return 
measurements from the D cup. N.B. the "Range*' is actually the 
heliocentric distance of the spacecraft 

In spite of the entry for 1982-103 there were photoelectrons 
observed in the El and E2 modes during the Uranus encounter 
during a spacecraft maneuver (day 24 of 1986). In addition, the 
threshold for detecting the solar wind in the EXT return A cup 
currents is apparently lower for Voyager 2 than for Voyager 1 
(refer to discussion above). As of day 160 of 1987, Voyager 2 
was 23.0 AU from the sun, and the first channel of the L mode 
(A cup) showed a DN of 132 in DC return. Channel 16 
registered a DN of 130 from photoelectrons. The "break" in the 
L mode spectrum occurs at the location of the solar wind peak 
in the A cup as indicated by the normal plasma data. 

To attempt to test the numerical loss calculations, we have 
concentrated on the Voyager 1 Cruise Science Maneuver from 
day 258 of 1978. From the PLS engineering output for this 
time of the maneuver, the electronics temperature (TPMQ 
12.6°C and the sensor temperature (TSEN) was -117.5°C. 
Modulation of the temperatures during the maneuver is clear 
from the engineering plots from this period (Figure 52). The 
times of the L mode spectra used, the DN value in channel 16, 
the corresponding current (corrected for temperature using 
interpolated values from Figure 9), and the angle between the D 
cup normal and the sun (using the nominal sensor orientation 
with respect to the spacecraft) are shown in the following table. 


I dc From CRSMVR 1978-238 

SCET 

DN 

I(fA) 

00, MM 

0132:25.276 

128 

< 10 s 

75.0° 

0134:01.276 

128 

< 10 s 

63.3° 

0135:37.276 

128 

< 10 s 

53.2° 

0137:13.276 

144 

2.50 x 10 3 

45.9° 

0138:49.276 

151 

3.25 x 10 5 

43.1° 

0140:25.276 

151 

3.25 x 10 3 

45.6° 

0142:01.276 

145 

2.58 x 10 3 

52.7° 

0143:37.276 

128 

< 10 3 

62.8° 

0145:13.276 

128 

< 10 3 

74.4° 


It is immediately obvious that the variation in current does not 
exactly track the variation in angle; whether this is related to the 
motion of the cup, shading of the instrument, bias in the timing 
information, hysteresis in the electronics or whatever is 
unknown. Another feature which is not understood is that when 
both the A and D cups are in shadow, at 0156:25.276 - angles 
between the sun and the normals of the A and D cups 66.6° and 
136°, respectively, negative currents were measured in both 
cups. At other angles, DN of 128 were recorded in shadow. 

To see if this behavior was peculiar to this maneuver (and to 
check for errors in the analysis), a cruise maneuver on Voyager 


ORIGINAL 
Of POOR 

2 was also examined in some detail. The first one, from day 
362 of 1977, was examined; the electronics temperature is not 
readily available so only the DN values from channel 16 (L 
mode, DC return, normal suppressor) are listed. The angles are 
obviously incorrect; the changes in the DN values are well- 
centered in the time period during which photoelectrons are 
observed in the El and E2 modes. Examination of the 
trajectory data indicates that the time labels are off by about the 
round-trip light time between the spacecraft and Earth (to 
within a minute). This error in the trajectory data precludes a 
better analysis. The trajectory data corresponding to the 
Voyager 2 Cruise Science Maneuver on day 155 of 1979 and 
the Voyager 1 maneuver on day 34 of 1979 are in Jupiter 
encounter format; the pointing information is supposed to be 
updated every 48 s, rather than every 4 m on the usual 
CRSMVR tapes (which is interpolated to provide the angles 
listed in the tables). 


loc From CRSMVR 1977-362 

SCET 

DN 

00, m 

0606:00.350 

128 

97.2° 

0607:36.350 

128 

109,0? 

0609:12.350 

128 

121.0° 

0610:48.350 

170 

130.0° 

0612:24.350 

177 

136.0° 

0614:00.350 

175 

137.0° 

0615:36.350 

169 

1310° 

0617:12.350 

140 

124.0° 

0618:48.350 

128 

113.0° 

0620:24.350 

128 

101.0° 


Examination of these maneuvers in more detail should provide 
better estimates on the angles through which the sunlight was 
entering the side sensor as a function of time. 


Predicted Photocurrent at the Collector Plate 

To compare the theoretical predictions of suppressor loss with 
the data, the quantity A^ of equation 12 must be computed for 
the illuminated area of the D cup collector plate. As photons 
are not refracted by the fields between the intervening grids as 
electrons are, Aloss is a one parameter function; the parameter 
can be taken to be the angle from the D cup normal at which 
sunlight enters the cup. To perform this calculation, the 
subroutine CPRESP was modified as shown by the handwritten 
notes in Appendix C. The actual code used is appended in 
Appendix E along with the effective loss area for both isotropic 
and cosine angular distributions of the secondaries. The optical 
transparency of the intervening grid wires is also listed. As the 
measured currents are from the positive ion mode, the lookup 
tables used are for -95.2 volts on the suppressor grid in the 
normal suppressor configuration. 

Comparison of Measured Currents with the Data 

The saturation photoelectron current density for gold under 
solar irradiation is j photo, am — 2.9 x 10~ 9 amperes cm' 2 [Grard, 
1973], For the Cruise Science Maneuver by Voyager 1 on day 
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258 of 1978 (the maneuver used by Barnett [1984] in testing the 
algorithm for the response function), this must be corrected for 
the heliocentric distance of the spacecraft, 4.152 AU. At 
Voyager, the expected photocurrent density from the illuminated 
part of the collector plate is jp/uu, a* cor = l- 7 x l®” 10 am P eres 
cm' 1 . Predicted currents as a function of angle, using this 
current density and values from the table in Appendix E are 
given in the following table; the loss current is in fA for both 
isotropic and cosine angular distributions. The illu min a t ed area 
is corrected for the grid transparency. 


loc Predicted for CRSMVR 1978-258 


SCET Bo.*- (lO 5 ^) (10 s fA) 


0135:37.276 53.2° 

0137:13.276 46.0° 

0138:49.276 43.3° 

0140:25.276 45.8° 

0142:01.276 52.9° 

0143:37.276 63.0° 


Note that the angles used are slightly different than those quoted 
above (pointing accuracy is no better than -0.02° in any case); 
the difference comes from a somewhat crude linear interpolation 
of the angles which are used in the table immediately above; the 
"better" values used in the table previously were not available 
when the predictive calculations were done. 

All of the predicted currents are larger than those actually 
measured; however, there is an additional effect which lowers 
the calculated predicted current. Sunlight on the suppressor grid 
will also produce photoelectrons; some of these (one can 
Actitnat* roughly half) will be accelerated toward the collector 
plate while the rest will be accelerated in the direction of the 
buckout grid and eventually lost to the sides of the cup (more 
details of where secondaries go after they are produced is 
discussed below). The grid wires are woven tungsten mesh. 
The expected photoelectron saturation current for tungsten at 1 
AU is bracketed by [Whipple, 1981] 

2.1 x KT* A cm' 2 < Vu*,, w < 8.1 x KT* A cm' 2 . Correcting 
for the distance of the spacecraft from the sun yields 

1.2 x KT 10 A cm' 2 < W* < 4.7 x KT 10 A cm' 2 . To have 

reached the suppressor grid, the light must pass through 6 grids, 
and to reach the collector plate from the suppressor grid another 
grid must be passed (which is at ground potential). The 
associated transparency can be estimated as 

(1 - l/[42 cos 43.3] ) 7 x (1 - ± ) 7 = 0.6693. The area of 

suppressor grid illuminated with sunlight is given by the overlap 
area of the aperture and the guard rings multiplied by 
1 - (1 - l/[42 cos 43.3] ) x (1 - -jj ) = 0.0558. The overlap 

area is 46.0 cm 2 ; hence, the area for emission is 2.56 cm 2 . 
Multiplying by the transparency yields an effective area of 1.72 
cm 2 , which is smaller than but not negligible compared to the 
effective area for photoelectron loss. 

If we take half of the average current density suggested for 


tungsten and combine it with this effective area, we obtain a 
negative current of 2.5 x 10 5 fA. 

From the table above, at 43.3°, the total predicted current is 
6 5 x 10 5 fA using the numbers for isotropic emission and 
2 JS x 10 5 fA assumin g the emission follows a cosine law. 
Given the uncertainty in the photoemission properties of the grid 
wires as well as the uncertainty in the fraction of photoelectrons 
emitted by the suppressor grid that will travel toward the 
collector, the predicted current (especially for the cosine 
distribution) is in excellent agreement with the maximum 
measured positive current of 3.25 x 10 5 fA during the Cruise 
Science Maneuver in question. 

We note that the DC return current should be negative in the 
reversed suppressor configuration due to the vastly increased 
suppressor efficiency. Unfortunately, no cruise full yaw turns 
were performed in the DC return mode with the suppressor 
position reversed; hence, this prediction cannot be verified. 

Other Sources of Photo- and Secondary 
Electron Currents 

Electrical Configuration of the Side Sensor 

Figure 53 is a scaled working drawing of the D cup showing the 
placement and support structures for die various grids. A 
comparably detailed drawing of one of the main sensor cups can 
be found in Barnett [1984] and in Barnett and Olbert [1986]. 
A less detailed drawing of the side sensor, in similar format, 
can be found in Barnett [1984]. The original version of that 
drawing is included as Figure 54. There are some differences in 
the dimensions which should be noted; they are probably too 
small to have any effect on any analysis done to date or 
contemplated for the future using the algorithm developed by 
Barnett [1984]. 

Although the dimensions in Figure 53 are quoted to the 
nearest thousandth of an inch (inches because those are the units 
in which the original mechanical specifications were made), the 
quoted uncertainties on all of the machine drawings for the 
various pieces of the assembly are typically ± 0.002 , if not 
larger. In labeling Figure 53, the intergrid spacings were 
computed from the nominal dimensions on the final machine 
drawings (including all of the Engineering Change Orders - 
ECOs) from the top surface of one grid to the top of the next 
(or to the top of the collector plate); the reference line is the top 
surface of the aperture grid. An expanded detail in the lower 
left comer illustrates the assembly of the 1 mil tungsten mesh 
(woven like nylon hose) to the grid support rings. There is a 
nominal 0.009” space between the top surface of the outer 
support ring and the top surface of the grid wires. This spacing 
has been included in calculating the between the grids and the 
top of the collector plate; it was not included in the dimensions 
on Figure 54. The mesh is stretched over each inner support 
ring and then secured with the outer ring (press fit). Excess 
mesh was then trimmed away at the bottom surface of both 
rings which were positioned to be flush. The nominal size of 
the openings in the mesh is 0.042 M x 0.042”. Variations from 
this are small and are documented in the appropriate engineering 
notebooks which record details of the grid assemblies. 

The thickness of the insulating spacers which electrically 


5.20 

3.37 

8.35 

4.85 

9.20 

5.19 

8.40 

4.88 

5.36 

3.45 

0.055 

0.047 


ORIGINAL PAGE IS 

OF POOR QUALITY 



- 11 - 


ORIGINAL PAGE IS 
OF POOR QUALITY 


isolate the bottom four grids were not specified exactly, rather, a 
range was given; to compute the spacing, the relative spacing of 
the centers of the holes for the anchor screws was used. These 
values were generally consistent with those inferred from the 
thicknesses of the grid support rings and the insulators. 

The intergrid spacings are listed in the table at the lower right 
in inches. The computed values are in the column labeled 
no min al. The column labeled Right 1 (i.e., SN002) are values 
measured after assembly of the cup using a microscope. The 
grid wires were brought into focus and the spacings computed. 
Most of the dimensions are fairly close and give a better idea of 
the range of error involved. The measured dimension F is 
almost certainly in error; not only does it differ from the 
computed value by 0.0049” (a lot), but also the difference 
between F and G is 0.0120”, which is smaller than the specified 
machine thickness of the outer support ring for the second 
suppressor grid. 

The intergrid spacings shown on Figures 53 and 54 are 
consistent to better than 0.002” for five of the spacings and not 
as consistent on the other three. The spacing E - D on Figure 
53 (based on the screw hole spacing) is 0.159” as compared to 
0.164” on Figure 54. This difference of 0.005” is annoying and 
probably small due to an error but follows directly from the 
dimensions to the anchor holes as given on the machine drawing 
for the lower housing. The spacing F - E is probably correct on 
Figure 53 and too large by 0.010” on Figure 54 (from including 
too much thickness for the shield ring). Finally, the spacing 
from the second suppressor grid and the top of the collector 
plate is 0.015" larger on Figure 54 than on Figure 53. I have 
triple and quadruple checked the distance from the aperture grid 
to the top of the collector plate; I believe that 2.338” is correct 
as shown on Figure 53. The collector plate thickness is 0.020” 
so this does not seem to be the source of the problem. The 
standoffs which separate the lower housing and collector plate 
come in two sizes 0.250” and 0.235”. There are two of each; 
the shorter ones also support solder lugs which come between 
them and the collector plate (these are used to connect the wires 
to the collector plate to make the electrical connection from the 
plate to the input amplifiers). If the shorter length were 
assumed for all of the standoffs, the grid to collector spacing 
would be larger by 0.015”; this is the likely source of the error 
on Figure 54. 

This latter discrepancy is probably the most important, as it 
would lead to more secondary electrons escaping from the 
collector plate than were found in the calculations discussed 
above. In other words, the suppressor is probably slightly more 
efficient than the numerical calculations done to date imply. 
This difference is probably not great enough to justify a 
recalculation of the tables for Au^ however. 

Overall, the aperture grid * to - collector plate spacing may be 
in error by up to 0.030” in the algorithms of Barnett [1984]. 
This maximum error is only 1.3% and certainly smaller than 
other uncertainties in the algorithm. 

Characteristics of the Buckout Grid 

The so-called "buckout" grid was included in the instrument 
design to fine tune the instrument to depress the electrical noise 
level. The effect of the buckout potential was not included in 


the algorithms of Barnett [1984] and Barnett and Olbert [1986]. 
Given the very small potentials involved and other uncertainties, 
such treatment is easily justified. However, the buckout 
potential may be important in understanding one of the features 
in the "photo spectrum”, and is not documented elsewhere in an 
easily accessible way. For reference then, we discuss it here. 

In the real instrument there is a certain amount of capacitative 
coupling between the modulator grid and collector plate, i.e., 
some field lines from the modulator grid "leak" through the 
grounded grids to the collector plate. As a result, the changing 
modulator voltage can produce a spurious current at the 
collector. The buckout grid is "tuned” to cancel this effect. 
The voltage applied to the buckout comes from the secondary 
winding of a transformer in the AC portion of the circuit for the 
modulator, hence, the buckout voltage is proportional to the AC 
part of the modulator voltage. The waveform is quasi-square 
and out of phase with the square wave voltage on the modulator 
grids. There are separate adjustments for the buckout potential 
for the main sensor and for the side sensor. The first failure (in 
1978) of the PLS instrument on Voyager 1 occurred in the 
buckout transformer for the main sensor, effectively degrading 
the sensitivity by a large amount which varied in time (the 
failure was "repaired” by thermally cycling the instrument, 
presumably cold soldering a connection which had come loose). 

After assembly, the voltage to the buckout grids was trimmed 
to minimize the modulator to collector feedthrough. For the 
16th channel of the L mode the buckout signal was measured to 
be 5 volts peak-to-peak for the main sensor and 2 volts peak- 
to-peak for the side sensor for SN001, the PTM or "Proof Test 
Model" which was flown on Voyager 2 (H.S. Bridge, private 
communication, July 14, 1987). These values can be scaled by 
the modulator voltage (AC component) to the other channels of 
the L mode and to the other operating modes (El, E2, and M) 
as well. 

It should be noted somewhere that SN003 was to be flown on 
Voyager 2 but was replaced with SN001 shortly before launch 
(the substitution was made June 17,1977 at ETR, i.e., the 
Eastern Test Range). Hence, SN003 was not flown on either 
spacecraft and is the "flight spare" in the group’s possession at 
M.I.T. 

A calculation of the capacitive coupling based upon planar 
grids rather than a mesh (and invoking other assumptions and 
approximations which are reasonable) is included as Appendix 
G. The solution follows the treatment presented on pp. 1235-7 
of Morse and Feshbach [1953]. 

Measured Currents in the Presence of Sunlight 
Angular Dependence 

If the angle between the D cup normal and the direction to the 
sun decrease to less than a critical angle when a negative 
potential is present on the negative modulator grid (refer to 
Figure 53), large modulated currents are measured at the 
collector plate in the El and E2 plasma modes. The spectra 
have a peculiar character in that some of the channels register 
threshold, i.e., 0 DN, regardless of the magnitude of the currents 
in the other channels. 

During 1977 and 1978 there were extended periods of time 
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during which sunlight entered the D cup at a fairly large angle, 
and these currents were measured in the electron modes for 
-100 days at both spacecraft during each year. The effect 
turned on and off suddenly in 1977 due to spacecraft rolls 
through relatively small angles; the cutoff occurred when both 
Voyagers were reoriented to acquire the DSN stations with the 
spacecraft high gain antenna (HGA). The turn on and turn off 
of the effect occurred due to very gradual changes in the D cup 
sun angle. As the spacecraft HGA stays pointed at the Earth, 
the angular separation of the sun and Earth on the sky change as 
the Earth orbits the sun; the amplitude of this variation 
decreases as the spacecraft recede from the Earth. In 1978 the 
photoelectron signal changed between the off and on states as 
the angle between the side sensor normal and the sun decreased 
from 78.5° to 77.5°. As the angle again increased through this 
range the effect shut off. This behavior was seen on both 
spacecraft, and is consistent with the range of angles over which 
the photo-induced currents have bee seen during the Cruise 
Science Maneuvers and special yaw maneuvers which have 
taken place as part of planetary encounter operations. It should 
also be noted that during 1978 while the effect was "on" the 
sun-D cup angle never decreased to less than 74°. In 1979, the 
minimum angle for the Voyager 2 spacecraft was 78.7° and no 
effect was seen; the minimum angle for Voyager 1 during 1979 
was larger. As both spacecraft have receded from the sun and 
the angular separation of the Earth and sun have decreased, the 
corresponding minimum angles have increased; the effect has 
not been observed again during cruise (all angles quoted are 
accurate to better than 0.5° and are probably accurate to 0.2° as 
found from observations of the steady solar wind during cruise 
maneuvers, H.S. Bridge, private communication, 1987). More 
details and documented evidence concerning the value of the 
critical angle can be found in Voyager Memorandum #156. 

Character of DC Return Spectra 

The DC return mode of the spacecraft has been referred to 
previously, and must be discussed again because they provide 
an important clue to the production of the observed modulated 
current Figures 55 and 56 show DC return spectra from the 
cruise phase for Voyager 1 and Voyager 2, respectively. 
During this phase of the mission sunlight was entering the D 
cup at less than the critical angle referred to previously. Both 
figures are a generalized version of the usual display of (the 
logarithm of) the current versus channel number (-logarithmic 
in energy per charge). The ordinate here is scaled to the 
logarithm of the magnitude of the measured current As the DC 
return currents can be positive or negative, regardless of the 
operating mode (i.c., the potential on the modulator grids), the 
plot format has been generalized to take into account both signs 
of the current. The abscissa is a similar generalization of the 
channel number. The basic unit is one channel width in the M 
mode. The width of an L or E2 channel is 8 basic units and the 
width of an El channel is 2 basic units (one E2 channel is 
spanned by four El channel widths where the two overlap in 
coverage). 

The "generalized data number 1 ' 1 through 8 corresponds to 
channel 16 of the E2 mode, 9 through 16 to channel 15 and so 
on, with decreasing magnitude of negative modulator potential. 


Data numbers 121 through 128 correspond to channel 1 of the 
E2 mode and numbers 127 and 128 correspond to channel 1 of 
the El mode; the "lower edge" of data number 1 is (a nominal) 
-5950 volts and the "upper edge" of data number 128 is -10 
volts. A 16 unit gap is included to represent the gap in 
coverage from -10 V to +10 V. Data number 145 then 
corresponds to channel 1 of the M mode (positive modulator 
potentials), and numbers 145 through 152 correspond to channel 
1 of the L mode. Channel 128 of the M mode corresponds to 
data number 272 and channel 16 of the L mode corresponds to 
data numbers 265 through 272. The "upper edge" of data 
number 272 corresponds to +5950 volts. 

The conversion between channel number and voltage is 
discussed in Bridge et al. [1977]; for the sake of completeness, 
we give it here. Let n be the channel number, which can take 
on values of 1 to 16 for the El, E2, and L modes and 1 to 128 
for the M mode. Let M be an integer with values of 8 for E2 
and L modes, 32 for El and 64 for M. Then the potential 0. at 
the lower edge ("bottom") of channel n, which is equal to the 
voltage at the upper edge ("top") of channel n + 1 (the upper 
edge potentials follow from considering "channel" 129 for the 
M mode and 17 for the other modes), is given by 

<►, = 60 10- /M - 50 (14) 

where 0„ is in volts. 

In Figures 55 and 56, data from the El, E2, and L modes are 
plotted; the times of the spectra are indicated across the bottom 
(the title "NEPTUNE MAGNETOSPHERE" is, of course, 
erroneous, and a glitch left over from rapidly converting some 
software to make these plots). Data from the M mode is not 
shown on these particular plots. The two traces for positive 
modulator potential are from the A and D cups, the D cup trace 
being at the positive threshold value. To plot the data as 
current, the conversion used is that of Figure 1 and Table 1. 
No temperature correction was made because it is not generally 
available. 

In the DC return mode, the PLS instrument works similarly to 
a Langmuir plasma probe but only in a superficial way. First 
because the aperture grid is grounded, the modulator potential 
can have no direct effect on the surrounding plasma; it can only 
effect plasma which drifts into the aperture. The exposed 
aperture grid is always at the spacecraft potential, so there is no 
way of using these measurements to directly infer the spacecraft 
potential by adjusting the modulator potential until zero current 
is measured. At large positive modulator potential, the positive 
ions should all be excluded, leaving a saturated negative current 
due to the electrons. At large negative potentials, all electrons 
should be excluded and a positive saturation current should be 
measured. This behavior does not usually occur in the DC 
return measurements although there are exceptions. 

For positive modulator potentials, the current saturates at a 
positive current which is due to the suppressor inefficiency 
(discussed previously). At lower positive modulator potentials, 
higher currents are measured as the solar wind is excluded less 
and less from the collector plate. For the times considered, the 
solar wind is more obvious on Figure 56 than on Figure 55 even 
though it is present in both. From the discussion in the 
previous paragraph, this curve should have started at a large 
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negative value, increasing toward zero as the positive solar wind 
is included with decreasing potential The curve is shifted up 
due to the suppressor potential In the solar wind, most of the 
negative current reaching the spacecraft is due to electrons with 
energies substantially less than the 95.2 volts (negative) on the 
suppressor; hence, in this case the suppressor acts as a rectifier 
allowing no negative current to reach the collector. There are a 
few cases, in the magnetospheres of Jupiter and Saturn when 
DC return currents were measured, that the electrons carrying 
most of the current through the aperture had energies in excess 
of 95 volts; in these cases, the expected behavior does occur; 
the total current to the collector plate is negative at large 
positive modulator potentials. 

For negative modulator potentials the expected behavior again 
does not occur. In this case, an increasing negative modulator 
potential should exclude more and more electrons while not 
affecting the ion current to the collector. At large negative 
potentials, all electrons should be excluded and the current 
measured should just be the positive saturation current due to 
the positive ions. In figures 55 and 56, increasing negative 
current is measured with increasing negative modulator voltage 
for the first channel of the El mode (data numbers 128 and 127) 
and the first 9 channels (data numbers 128 through 56) of the 
E2 mode. Due the way that currents are synchronously detected 
in the normal plasma mode [Bridge et al. t 1977], the currents in 
the first 9 channels of a regular plasma measurement in the E2 
mode (and the first channel of El) will not be detected because 
the current is out of phase with the detector, a level of 0 DN, 
i.e., threshold, will be entered into the telemetry stream. The 
modulator potential at the upper edge of channel 9 is -750 volts. 
This corresponds to an energy close to the maximum in the 
secondary emission yield curve for gold, an important clue to 
the origin of this effect At higher negative potentials, the 
negative current is decreasing in magnitude, but is still large and 
negative at -5950 volts, a potential sufficiently large to easily 
exclude all ambient electrons. 

The switch in suppressor potential between -95.2 volts in E2 
and -8.05 volts in El introduces a subtle change. The two traces 
from El and E2 are different just below -95 volts (data number 
less than 103), although they are approximately continuous at 
the "top” of the El mode. Another effect is present also in that 
channel 3 of E2 (data numbers 111 to 104) corresponds to 
nominal potentials from -56.7 to -92.3 volts. However, 
electrons in this energy range cannot pass the -95.2 volt 
potential barrier of the suppressor grid, so a current with the 
modulator at this channel is also a mystery. These changes 
imply that both the modulator and suppressor potentials 
determine the current measured at the collector plate. 

The position of the suppressor potential also affects the 
currents. For both Voyager 1 and Voyager 2 there are sets of 
EXT return spectra acquited within two hours of those spectra 
shown in Figures 55 and 56. Figures 57 and 58 show DC return 
spectra for Voyagers 1 and 2, respectively, in the reversed 
suppressor configuration. In both cases, there are qualitative 
changes in the shapes of the curves and the current is larger in 
magnitude by roughly 50% in the reversed configuration. This 
behavior suggests that in the normal configuration, the 
suppressor excludes more current from reaching the collector 


plate. 

In Figures 59 through 69 are EXT return scans from the 
Voyager 1 encounters with Jupiter and Saturn. They illustrate 
(sometimes) the "correct" behavior of the spectra and also some 
of the difficulties of interpretation. Figure 59 is from the 
middle magnetosphere of Jupiter. Photoelectron current 
dominates the A cup signal. E2 and D are continuous (both use 
the same suppressor voltage); note that E2 exhibits an increasing 
(positive) current with decreasing (negative) potential in the 
lowest 4 channels. This out of phase behavior is similar to that 
produced by photoelectrons and is, presumably, due to 
secondaries produced by hot energetic electrons in the ambient 
medium. The side sensor was in not illuminated with solar UV 
at this time (UV reflected from Jupiter might be able to enter 
the cup - I haven’t checked the geometry - but is probably of 
too low an intensity to have an effect). Note that the transition 
through zero net current occurs in the El scan. That it does not 
occur at zero modulator potential is not surprising; the 
spacecraft is not illuminated uniformly with the positive ion and 
electron fluxes, and this will be manifested in the zero current 
potential. The "crossover’' is in no way related to the floating 
potential of the spacecraft with respect to the plasma due to the 
grounded aperture grid. 

Figures 60 through 64 contain data from the warm Io torus. 
Effects of photocurrent loss in the A cup, hot electrons, and the 
change in the ram direction of the positive ions are all present, 
as well as what appears to be the effect of rapid time variations 
resulting in time aliasing of the A cup data in Figure 61. 

Figure 65 is illustrative of (but not the best example of) data 
from the cold torus, and Figure 66 is from the torus outbound 
just prior to the encounter with the Io flux tube. Figure 67 is 
the from the period of solar occupation; at large positive 
modulator potentials the A and D cups track each other well as 
the shadow of the planet has cut off the photoelectron 
production in the A cup. With both the A and D cups looking 
away from the positive ion flow direction, the positive ion flux 
reaching the spacecraft is poorly sampled. As a result the 
transition through zero current occurs at a relatively large 
negative modulator potential (the entire trace should be shifted 
up) and the A and D cup spectra are relatively flat (neither cup 
is pointed in the correct direction to sample ions which a 
positive modulator potential would exclude). Figure 68 shows 
the first set of spectra back in sunlight; the A cup is still low, 
presumably due to a large electron flux. 

Finally, Figure 69 is one of the more interesting ones from 
the Saturn encounter. It is from the outbound crossing of the 
Dione - Tethys plasma torus. The positive potential traces 
appear to have the wrong sign for the current; at least that 
would make more sense than the data as displayed. The effect 
of the suppressor voltage in excluding electrons is illustrated 
nicely by the El and E2 data in this case. 

In principle, it should be possible to use plasma parameter 
determinations from the normal plasma data in conjunction with 
the DC return spectra along with the condition of zero current at 
zero modulator potential (equivalent to modulator off EXT 
measurements, of which we have none in the magnetospheres) 
to estimate the net current to the spacecraft due to the 
(subsonic) electrons and ions with energies per charge above 
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5950 volti. In practice this would be difficult because the 
electron plasma and DC measurements would have to be 
collected for contamination by secondaries within the housing 
itself. More troublesome than that is the fact that the secondary 
electron emission characteristics of the materials from which the 
spacecraft is made are not really known (discussed above in the 
section on photoelectron yields). 

Electric Fields Within the 
Side Sensor Housino 

For the purpose of calculating the instrument response to a 
variety of plasma conditions, Barnett [1984] and Barnett and 
Olbert [1986] make the simplifying assumption of uniform 
electric fields between the various grid planes. For most 
purposes this is adequate. Further refinements, such as ta kin g 
into account the nonuniformity of the field in the immediate 
vicinity of the grid wires, effective makes the problem 
enormously more complicated, if not impossible due to the 
small differences between the sensors as designed and as 
actually configured once assembled. An example (the subject of 
grid spacing) has been mentioned above. There are two 
examples where the effects of the field nonunifonnities may 
play important roles in the effects of photoelectrons and 
secondary electrons within the instrument housing. We can also 
give some quantitative estimates of these nonunifonnities, if we 
make some approximations. The regions in question are in the 
vicinity of the modulator grid support ring and in the vicinity of 
the suppressor grids. 

Fringing Electric Field Near the Negative Modulator Grid 
Support 

In light of the evidence presented, the currents produced when 
the PLS side sensor is appropriately illuminated with sunlight 
must be due to photoelectrons produced on the inner surface of 
the inner support ring of the negative high voltage modulator 
grid. Referring to Figure 53, it does not appear possible for 
such electrons to reach the collector directly due to the presence 
of the shield ring. A more important objection to a direct 
collection of these electrons is the fact that all electrons 
produced at this location will be "born” at the modulator 
potential; photoelectrons are cold, with an average energy of 
less than-1 eV, small compared to the other potentials which 
are present To assess the possibility of a direct contribution to 
the current, it is necessary to consider the electric field 
configuration in the vicinity of the ring. 

The actual boundary conditions for the appropriate Dirichlet 
problem are rather complicated even though there is cylindrical 
symmetry. To simplify the problem so that it is more easily 
treated numerically, we can replace the true boundary value 
problem with a two-dimensional analogue. If the scale length 
of changes in the potential is small compared to the support ring 
radius, this approximation should be adequate for our purposes. 
Figure 70 shows the boundary conditions for the two- 
dimensional problem. The electrostatic potential is sought for 
the indicated surface at some potential V^; all other surfaces 
are at ground potential. Equipotentiais were calculated 
numerically and are displayed in Figure 71, and electric field 


lines in the near vicinity of the support are displayed in Figure 
72 (from a Laplacian problem solver available on Project 
Athena, C. K. Clark, private communication, June, 12, 1987). 
The modeling used a fairly coarse grid spacing but does indicate 
that the "fringing" effect around the support rings dies out after 
a displacement toward the center of the cup equal to about twice 
the width of the combined inner and outer ring supports. This 
is, also roughly, half of the distance from the inner surface of 
the inner ring to the inner edge of the shield ring. 

An approximate solution to the problem is available in 
analytic form using conformal mapping techniques. This 
approximation includes only the inner wall of the support ring 
with their surface facing the shield ring extended to infinity. 
The appropriate map (inverted) is from Kober [1957] and is 
shown in Figure 73. 

A detailed calculation of electron trajectories from the surface 
of the support to the the plane of the shield ring has not been 
carried out; however, it appears unlikely that the fringe field can 
impart enough of a "kick" to the electrons to enable them to 
clear the shield ring. A proper calculation would need to take 
into account the drift space between the various ground grids 
and buckout grid. The intergrid spacing is small enough and the 
modulator frequency low enough that the problem can be treated 
as a static one with photoelectrons first produced at the upper 
potential and then at the lower potential for a given modulator 
step. 

The problem is simplified somewhat by the low energy of the 
photoelectrons. As the energy the electrons gain from the 
modulator potential is large compared to the energy with which 
they are "bom” and the buckout potential is a small fraction of 
the modulator potential, the trajectories of the photoelectrons are 
independent of the modulator potential. This statement is true 
only above die ground grid above the first suppressor grid in the 
normal suppressor configuration. As discussed immediately 
below, the suppressor potential, which is independent of the 
modulator potential, produces a nonuniform field in this region. 
The effect of this constant field will tend to deflect toward the 
cup wall the electrons coming from the modulator grid. The 
deflection will then determine where on the shield ring the 
electrons impact; the deflection should decrease with increasing 
modulator potential. 

Further Work on the Fringing Problem 

After the initial assessment of the "fringing field" problem 
discussed in the section immediately above, it became obvious 
that the question of whether electrons produced on the relatively 
large surface area of the modulator support ring can indeed 
reach the collector plate. As of January, 1988, an experimental 
answer based on the characteristics of the Flight Spare PLS 
instrument is still unavailable. In December, 1987 an extensive 
set of calculations was undertaken (by G. S. Gordon, Jr.) to find 
a theoretical answer to this question. Details and results of this 
calculation are outlined in this section (note that the figures 
referred to herein are out of numerical older with some of those 
referred to in following sections). 

The problem is solved in a two-dimensional approximation 
using the map from Korber [1957] depicted in Figure 73 (this 
mapping function is also given in problem 10.21 on p. 1315 of 
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Morse and Feshbach [1953]. The appropriate Schwarz- 
Christoffel transformation is specified by (note that the 
numbering of the equations is out of sequence) 


dz _ k'Jw-l 
dw n wVw - c 


( 22 ) 


with k a he 1/2 and c > 1, k > 0. The solution is 


z (w) = — cosh 
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2w - c - 1 


c - 1 
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n sc 


[c + 1] w - 2c 


[c - 1] w 


(23) 


In this problem Jk = 0.452 inches, h = 0.312 inches, and 
c- 2.0988; the dimensions again correspond to the nominal 
values given in Figures 53 and 75. Details of the computer 
routines used to make the plots discussed immediately below 
and to obtain w (z) from eqn. (23) can be found in Appendix F. 

Contour maps of Re [z (w)] are shown in Figures 78 through 

81 for different regions with contour lines at 0.1 and 0.01 
levels. Corresponding maps of Im [z (w)] are shown in Figures 

82 through 85. The inverse mapping w (z) was found using a 
variation on the Newton-Rap son technique (see Appendix F and 
the appropriate source listings for further details). Contour plots 
of Re [w (z)] and Im [w (z)] are shown in Figures 86 and 87, 
respectively. The contour levels are powers of 2 in these two 
figures. In Figure 86, heavy lines indicate the contour levels of 
-1, 0, and +1; heavy lines in Figure 87 show the contours at 0 
(along the periphery) and +1. Note that the rectangular region 
in the upper left hand corner indicates the modulator grid 
support. Dimensions are given in inches in these figures with 
the 0.0 along the vertical scale located at the ground grid just 
below the negative modulator grid (cf. Figure 53). 

Given the mapping function, Laplace’s equation can be 
solved in the w plane and the solution mapped into the z plane. 
From p. 1245 of Morse and Feshbach [1953], we choose the 
complex potential 

y 

F a y + i x = t- In w (24) 


which gives zero potential on the positive real w axis and 
potential V on the negative real w axis. Hence, 

w(F) = e* f,v (25) 

and contours of constant Re F = y are lines of constant 
potential which can be mapped to the z plane. These are 
shown in Figure 88. Contours of constant Im F ~ % give the 
electric field lines; these are shown in Figure 89. 

Given the electrostatic potential of Figure 88, we assume 
photo- or secondary electrons are "bom" on the inner edge of 
the modulator support ring. For a given set of initial conditions, 
the trajectories of these electrons can then be integrated to the 
plane containing the shield ring to see what fraction pass the 
ring and can continue to the collector plate (this, of course, 
assumes that the electron energy is greater than that expended in 
climbing the potential barrier of the suppressor grid. To make a 
set of representative calculations, we have assumed that the 


suppressor potential is zero (representative of the reversed 
suppressor configuration) and that the electrons are emitted 
normal to the surface. Figures 90 and 91 show the resulting 
trajectories for -10 volts on the modulator grid and initial 
electron energies of 0 eV and 2 eV, respectively. Figures 92 
and 93 show trajectories for modulator potentials of -6000 volts 
and the same two initial electron energies. 

On Figures 90 through 93, the inner edge of the support ring 
is in the upper left comer. Small tic marks on the y axis locate 
the two ground grids and buckout grid located at dimensions C, 
D, and E on Figure 53. For the purposes of this calculation, the 
buckout grid was also assumed to be at ground potential. As 
the electrons have been accelerated through the modulator 
potential before reaching it, this potential should produce a 
negligible effect on the trajectories. The edge of the shield ring 
is indicated by the rectangle on the lower left Again all 
dimensions are in inches. The information of the previous 
figures is summarized in Figure 94. The equipotential contours 
in the region between the negative modulator grid and the first 
grounded grid are shown at the top; the trajectories shown in 
Figures 90 through 93 are shown below with the physical region 
covered to the same scale. 

The immediate and obvious conclusion to be drawn from the 
figures is that the electrons get enough of a sideways ’ kick" to 
make it past the shield ring. Therefore, such electrons will reach 
the collector plate and contribute a negative current. This 
current is a constant, i.e., has no modulated component To 
attempt to determine the origin of the modulated component 
detected in the normal plasma mode, we must now consider the 
effects of the suppressor potential. 

Electric Fields in the Space Above the Suppressor Grid 

In both the reversed suppressor configuration and in the 
normal suppressor configuration, the field line configuration 
within the cup is somewhat complicated due to the details of the 
cup design. In the usual analysis algorithms, the electric field 
lines are all assumed to be parallel with the symmetry axis; this 
approximation is not necessarily as good in dealing with the 
effect of electrons produced on the modulator support ring. In 
the reversed suppressor configuration, the previously computed 
suppressor efficiencies will be modified by the fact that edge of 
the collector does not extend as far as the edge of the second 
suppressor grid itself (see Figure 53); there are modifications in 
the vicinity of the collector plate edge because, past the edge, 
the suppressor potential must drop across a greater distance, i.e., 
the distance to the bottom of the cup housing which is also at 
ground potential. Between the second suppressor grid and the 
first suppressor grid (clamped to ground in the reversed 
suppressor configuration) the field lines should indeed be 
straight (ignoring all of the incalculable second order effects 
present in the real cup). It is important to note that the 
(approximate) plane occupied by the shield ring and the first 
suppressor grid is at ground potential in this case. 

In the normal suppressor configuration, the second suppressor 
grid is held at ground potential, and the field lines between the 
second suppressor and first suppressor grid are straight, to first 
order. A complication in the field structure now arises in the 
space between the plane occupied by the first suppressor and 
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shield ring on the one hand and the ground grid which separates 
the first suppressor and buckout grids on the other. 

The space in question can be modeled as a circular cylinder 
whose sides and one end are at a potential of zero and whose 
other end is at a finite potential out to some radius a (the hole 
in the shield ring) and held to zero potential from radius a out 
to the edge of the cylinder at radius b (the space above the 
shield ring). The configuration exhibits sufficient symmetry that 
the potential can be written explicitly as a Fourier- Bessel series. 
Choose z = 0 as the plane of the ground grid just below the 
buckout grid, and let positive z be in the direction toward the 
collector plate; the boundary condition of = 0 at z = 0 is 
satisfied by taking the potential to be proportional to sinh kz. 
The combination of cylindrical symmetry and finite potential at 
r = 0 imposes a radial dependence of J Q {kz). The boundary 
condition <X> = Oatr=bis satisfied by taking k = x^ , lb where 
„ is the n th zero of the Bessel function of order zero. Let the 
spacing between the plane of the ground grid and that of the 
shield ring be denoted by L. Again, we ignore the 0.024" 
difference between the top surface of the shield ring (gold plate 
on a copper-clad piece of epoxy glass) and the top surface of 
the first suppressor grid. Then using the expansion 

*) = £ A m J 0 (ks) sinh kj 

* ■ l 

where k. m x^ t , fb and applying the boundary condition 
<t> (r, L) = V for 0 < r < a 

and zero otherwise, we can find A m . After some algebra 

. _ „ „ a csch kj, ■/, (k^a) 

A - m2 

The series does not converge very rapidly (for n * 1 to 4, A m 
is +10.7, -0.53, -1.0, and +0.99, respectively). 

The same calculation can be made for a two-dimensional 
geometry, i.e., eqn. (15) is now replaced by 

4> (r, z) = 2 A * cos (M sinh V (15a) 

*add * 1 

with k* m rm /2b. Applying the boundary condition eqn. (16) 
for this formulation yields 

4V csch 

A m a si n(k*r) (17a) 

it n 

which is an equally unenlightening result. 

An integral formulation can be found by letting b -> An 
analysis of this case yields 

4* (r, z) = V | J 0 (kr) J } (ka) d (ka) (18) 

which is more elegant but no more illuminating. This case can 
also be approximated by a two-dimensional analogue. In this 
case we obtain 


(15) 

(16) 
(17) 


<l>(r,!) = 2v] c <M*x»infc» (19) 

rc i k sinh kL 

This integral can be evaluated in terms of simple functions by 
noting that cos kx sin ka = l /2(sin k (x + a) — sin k (x - a) and 
using 4.114.1 of Gradshteyn and Ryzhik [1980], viz. 


r sin a k 

*** dk ~ tan' 1 

♦ n z 

tan 

* i* a 

tanh — - 

ft * 

sinh kL 

2 L 

2 L 


A plot of the field lines is shown in Figure 76 (note that the 
numbering is off. Figures 74 and 75 have not yet been 
introduced). The fact that the field lines are not really parallel to 
the symmetry axis of the cup is not included in either the 
algorithm of Barnett [1984] for the side sensor or in the 
numerical integration scheme PMODEL. Note that the effect on 
charged particles entering the cup will be to focus ions and 
defocus electrons. Put differently, the effect is to make the hole 
in the shield ring larger for ions and smaller for electrons with 
the effect becoming less important with increasing particle 
energy. 

To better assess the potential configuration, eqns. (15) and 
(17) were numerically evaluated by G. S. Gordon, Jr. to produce 
a contour plot of $ (r, z). This is shown in Figure 77; note that 
the figure is inverted with the suppressor grid at the top and the 
ground plane at the bottom. The dimensions are in inches, and 
the plot space (0.165 x 2.765) corresponds to the dimensions on 
Figure 75 (the vertical dimension used in Figure 76 does not 
include the 0.015 inch thickness of the shield ring). The 
contours are at 10% intervals, except for those nearest the 
boundaries, which are at 1% above/below the boundary value. 
The mesh used is 999 x 999 and the sum of eqn. (15) was 
terminated at either the 999th term or when the next term gave 
a further correction of less than 1CT 6 , whichever came first 

At this time, a rigorous, quantitative analysis of the difference 
between switching between normal and reversed configurations 
cannot be given. Qualitatively, one might think that the 
reversed configuration should lead to more electrons in the 
vicinity of the hole in the shield ring actually passing through 
the hole. This is consistent with the increased negative current 
observed in the EXT return mode, reversed suppressor 
configuration, when sunlight is illuminating the negative 
modulator support ring. 

From early in the mission of Voyager 2 on day 243 of 1977 
(11 days after launch), there are DC return spectra for both the 
normal and reversed suppressor configurations. At this time the 
angle between the sun and the D cup normal was 73.3°, hence, 
the negative modulator support was illuminated but not the 
shield ring (and suppressor grid) or collector plate. A plasma 
measurement from 1977-238:2011, normal suppressor, 
modulator on, shows the "turn on" in channel 10, i.e., channel 9 
is at threshold, as are channels 3, 4, 6-9, 15, and 16. This is in 
sharp contrast to the turn on location in the PLOMAN 2 spectra 
(see below). A corresponding DC return measurement on 
1977-249:1340 indicates a maximum (negative) current of 
-7.04 x 10 6 fA in channels 9 and 10, i.e., a maximum at -750 
volts (negative) on the modulator grid. This is consistent with 
the turn on in one of the "voids" in the Jovian magnetosphere 
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(see below). An earlier DC measurement at 1977-243:1856, 
gave a broad maximum of -6.40 x 10 6 fA in channels 7 through 
11, indicating that the maximum is flat to -10% (resolution in 
this current range for the DC currents, cf. Table 1). 

About an hour after the DC measurement in the normal 
configuration, at 1977-243:1954, there is a DC return 
measurement in the reversed suppressor configuration. Here the 
maximum current is —9.60 x 10® fA and extends from channel 6 
through channel 10. This is an increase of 50% in the current; 
the fractional increase in other channels is less, ranging down to 
% for channel 16. The only logical explanation is that 
increased refraction of secondary electron trajectories in the 
normal configuration prevents more electrons from reaching the 
collector. Given that additional refraction of electrons emitted 
from the modulator support ring should direct them farther 
toward the center of the cup, this is also consistent with most of 
the electrons originating on the underside of the modulator grid 
itself, presumably due to internally reflected photons. 

It is interesting to note that the next available plasma 
measurement in the reversed suppressor configuration, at 1977- 
242:2332, shows a turn on in channel 8, i.e., at a modulator 
voltage of 400 volts. This is consistent with the movement 
downward (in modulator potential) of the broad maximum in the 
DC return current; it is not simply consistent with the yield 
curve for any exposed materials in the interior of the cup. We 
return to this point below in considering the fractional 
modulation in the presence of sunlight. 

A Quautative Explanation of the Effect 

Illuminated Area of the Support Ring for the Modulator Grid 

To quantitatively evaluate the spectra due to photoelectrons, we 
need to know the illuminate area of the inner surface of the 
inner ring supporting the negative high voltage modulator grid 
as a function of the angle of illumination. The geometry is 
illustrated in Figure 74. Let the inner radius of the ring be a (- 
2.457” from Figure 53) and let the height of the ring be h (- 
0.140"). Assume the ring is tilted at an angle 0 with respect to 
the viewing direction. The illuminated area (in projection) 
appears as the intersection of two ellipses of semi-major axis a 
and semi-minor axis a cos 9 whose centers are displaced by an 
amount h sin 0. The geometry is illustrated in Figure 74. Note 
that the illuminated area is conveniently divided into Regions A 
and B as indicated in the figure. The point where the lower 
branch of the upper ellipse y x (x) coincides with the upper 
branch of the lower ellipse y % (x) separates the two regions; its x 
coordinate is given by xq = — ( h tan 0 fl). 

A straightforward, but tedious, integration gives the projection 
of the illuminated area Am**. Define the parameter 
P s h tan 0 /(2a). The illuminated area is then given by 

Am** = tih sin 0 Vl - P 2 + 2 a 1 cos 0 sin" 1 P (21) 

For small 0, A**, -► 2ah sin0 as expected. 

We have investigated the variation with angle of die 
modulated part of the current in the E2 spectra acquired during 
the PLOMAN 2 maneuver within the Jovian magnetosphere. 
The following table lists the times of the spectra, the angle 0, 
the illuminated area of the modulator support from equation 21, 


and the peak current For the first two spectra the peak current 
appeared in channels 10 and 11 (i.e., from 750 to 1373 volts); 
the peak current was located in channel 11 for the rest of the 
spectra (1017 to 1373 volts). 


E2 Spectra From PLOMAN 2 Maneuver 1979-062 

SCET 

6 

AtUMN cm 

u (10 4 fA) 

1825:19.471 

75.3° 

4.28 

1.59 

1826:55.471 

63.6° 

3.97 

2.63 

1828:31.471 

53.4° 

3.56 

3.15 

1830:07.471 

46,1° 

3.19 

3.26 

1831:43.471 

43.2° 

3.04 

3.15 

The change 

in illuminated 

area of the 

support ring is 


obviously not the major factor in causing the current to increase 
as 0 decreases; the illuminated area is actually decreasing due to 
projection effects. If we assume that the photoelectron current 
varies as the secant of the angle between the incident photons 
and the normal to the support ring surface, we still cannot 
account for the factor of two increase in current as the cup is 
rotated more into the sunlight 

Let a be the angle of incidence of the photons with respect to 
the normal to the inner surface of the ring. Then 

< Aiu ZL> = 2ah (— - sin' 1 p) + 4a 2 cot 0 (1 - Vl - P J X26) 
cos a 2 

where p is defined as above and < * * > indicates an average 
over the illuminated angle. In this case we obtain 


E2 Spectra 

From PLOMAN 2 Maneuver 1979-062 

SCET 

e 

AiUwm _ 2 

cos a 

<10 4 fA) 

1825:19.471 

75.3° 

6.73 

1.59 

1 826:55.471 

63.6° 

6.84 

2.63 

1828:31.471 

53.4° 

6.89 

3.15 

1830:07.471 

46.1° 

6.91 

3.26 

1831:43.471 

43.2° 

6.91 

3.15 


The averaging over cos a increases the effective illuminated 
area by about 50%; dividing A^ by sin 0 gives a roughly 
constant value of 4.4 cm" 2 . These variations show that even 
including the angular dependence of the electron yield per 
photon, the variation cannot be accounted for. 

To produce a modulated current, electrons produced at the 
modulator potential must be involved. A final source for such 
electrons is the modulator grid itself illuminated by photons 
reflected from the inside of the cup (E. C. Sittler, Jr, private 
communication, December, 1987). Photoelectrons produced on 
the underside of the grid wires will be accelerated toward the 
collector plate while those produced by direct illumination on 
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the top will tend to be accelerated toward the aperture grid and 
leave the cup. 

In the review on reflectivity in the ultraviolet by Koiler 
[1965] p the reflectivity of thin gold film is -20% for photons of 
energies between 6 and 10 eV. Such photons form the bulk of 
the solar spectrum in the energy range for which photoelectron 
yields are appreciable for many materials [ Grard , 1973]. 
Although photoelectron yield versus photon energy is not 
readily available for tungsten, it is reasonable to assume that 
-1/5 of the photons incident on the inner surfaces of the cup 
will be reflected back to illuminate the negative modulator grid 
wires. In this case the illuminating flux goes like cos 0. At 
normal incidence, the grid wire area is 
na 2 [1 - (1 - ) 2 ] = 5.76 cm -2 , i.e., comparable to the directly 

illuminated area of the modulator support ring. At angles of 
incidence of 75.3° and 43.2°, respectively, the areas of the 
modulator grid wires illuminated by sunlight are 1.46 and 4.20 
cm* 2 , respectively. These values should be reduced by a factor 
of 0.2 (the reflectivity coefficient). Adding the values from the 
table apparently again give too low a modulation with angle. 

One final correction is the transparency of the aperture grid, 
positive modulator grid (grounded), and negative modulator grid 
to the photons hitting the support ring. The transmission 
probability can be estimated as 

r=(1 -~k scc9)1 < 27 > 

This effect also yields too small a modulation with angle of 
incidence: in changing the angle of incidence from 75.3° to 
43.2°, respectively, the transparency increases from 0.692 to 
0.842 if estimated from eqn. (27). 

It is important to note that as the PLOMAN 2 maneuver 
occurred within the magnetosphere of Jupiter, the positive ions 
were flowing about at right angles to the spacecraft- sun line 
[Sands and McNutt , 1988]. This is not an important 
consideration for the plasma mode measurements as such; 
however, it identifies the problem with attempting to determine 
the angular modulation of the unmodulated electron current 
produced by photoelectrons. This component is only available 
from DC return measurements during a period in which the 
spacecraft is rotating. Such measurements are made during the 
Cruise Science Maneuvers mentioned earlier. These maneuvers 
are made only in the solar wind, so the angular modulation of 
the photocurrent also includes the positive contribution of the 
current due to the solar wind at the smaller angles of incidence 
(<-63° at which the solar wind can reach the collector plate - 
cf. Figure 53). In addition, there is the positive current which 
results from photons on the collector plate and escape due to the 
inefficiency of the suppressor grid, an effect discussed in detail 
above. We can attempt to disentangle these effects by 
differencing DC return currents from the L and E2 modes in the 
same channel acquired at about the same angle (there is some 
difference in angle due to the finite rotational speed of the 
spacecraft and the non-simultaneity of the L and E2 
measurements). 

We again consider the Cruise Science Maneuver on Voyager 
1 from day 258 of 1978. If we use the nominal DN to current 
conversion of Table 1, the currents are smaller than those listed 


in the above table in the text. It does appear that the time delay 
between L and E2 acquisition is too slow compared to the 
angular change in spacecraft orientation to use this technique to 
come to a definitive conclusion about the angular variation of 
the DC photocurrent. 

The variations are plotted as a function of angle in Figure 95. 
The modulated current is fairly consistent with a cosine 
variation, except for the turnover at angles of incidence of less 
than —50°. The DC curve indicates an even larger variation at 
larger angles, but this could be due to not correctly 
incorporating the electronics temperature in the DN to current 
conversion. The obvious conclusion to be drawn is that the 
internally reflected photons illuminating the underside of the 
negative modulator grid must play a major role in producing the 
photoelectrons in the instrument. In light of the lack of 
information in the literature on the photoelection yield of 
tungsten as a function of photon energy, as well as the 
geometric complexities of internal reflection within the 
instrument, it is likely that a quantitative model of the 
photoelectron response could only be derived from extensive 
laboratory calibrations of the flight spare (under appropriate 
conditions). If back illumination of the modulator grid is the 
major source of photoelectrons, then their utility as tracers of 
secondary electrons produced inside the instrument must also be 
considered limited. 

Production cf Intermediate Secondaries 

Given the photoelections produced at the support ring and/or 
on the modulator grid itself, we turn to the question of how this 
can give rise to a modulated current at the collector plate. As 
the modulator potential increases, the photoelectron current 
produced will increase somewhat but reaches a saturation value 
at an applied potential of a few volts. This could explain the 
initial increase in current with negative modulator potential, but 
cannot explain the peak at -750 volts. The location of the peak 
current in the DC return spectra during the Cruise Science 
Maneuver on day 258 of 1978 is located in channels 8 and 9. 
PLOMAN 2 spectra shows a "turn on" of the effect in channel 
9, which is consistent with the location of the peak in the DC 
return spectra. The thresholds of channels 9 and 10 are 550 and 
750 volts, respectively, so 750 volts is certainly past the 
maximum. According to Kollath [1956], the peak in the 
secondary yield curve for gold is 800 volts and for tungsten 650 
volts. Hence, the location of the peak in current versus 
modulator potential is consistent with secondaries produced 
from a tungsten surface, i.e., from the grid wires, but not from a 
gold surface, i.e., the plating on the cup surface. 

An exception to this occurs on Voyager 2 during its 
encounter with Jupiter in 1979 on day 190 from -1150 to 
-1155. This period is during the crossing of one of the "voids" 
identified in the Voyager plasma data [McNutt et al. t 1987; 
Khurana et al. t 1987]. During this period, the E2 current is 
apparently out of phase with the synch detector for all channels 
through channel 9, i.e., the "turn on" occurs in channel 10 
above its threshold of 750 volts. This suggests that at very high 
electron temperatures, the secondary electrons produced off of 
the gold plating can produce a significant effect 

Secondary electrons capable of reaching the collector plate 
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must have an energy in excess of -95 volts to pass the potential 
barrier established by the suppressor potential (secondary 
electrons can also be produced on the grounded grid between 
the suppressor and collector in the normal suppressor mode). 
This implies that the electrons must originate on either the 
modulator grid itself (or the support ring) on on the suppressor 
grid. The fact that sunlight produces a modulated current means 
that this current must be produced on a surface whose potential 
varies with respect to the modulator as the modulator voltage is 
stepped. By elimination, this surface must be the suppressor 
grid itself. At the peak current, the suppressor is at a potential 
of 0 - 95 * -650 volts with respect to the modulator; this, of 
course, corresponds to the energy of the peak in secondary yield 
from tungsten. 

Most of the secondaries produced at the suppressor grid will 
be produced on the sunlit, or outward facing, side (along with 
photoelectrons which constitute another unmodulated 
component). These will be accelerated by the suppressor 
potential out of the cup; however, if the modulator voltage is 
greater than the suppressor voltage (as is the case except for the 
lowest 3 channels of the E2 mode), these will be reflected back 
into the cup by the modulator potential. Furthermore, they will 
have just enough energy to climb over the potential barrier 
between the suppressor grid wires, which is fractionally smaller 
than that on the wires themselves. Only those electrons initially 
emitted nearly parallel to the cup axis will not drift to the sides 
of the cup and be lost; however, those that do make it through 
will reach the collector plate. For a given lighting condition, 
the number of electrons striking the collector will be 
proportional to the secondary electron yield of tungsten at an 
energy corresponding to the difference between the modulator 
and suppressor voltages. Some of the possible electron 
trajectories and secondary production sites are illustrated 
qualitatively in Figure 96. 

The other case of real interest concerns the production of 
secondary electrons within the instrument by energetic 
secondary elections in the local plasma environment Again, to 
produce a modulated signal (which can either add to or subtract 
from the modulated primary electron current) the modulator 
potential must be involved. As noted by Vasyliunas [1968], a 
component of secondary flux in phase with the primary flux 
(and therefor simply adding to it) will be produced at both the 
suppressor grid (proportional to the yield at energy 
Ephmar, ~ eV ,q>) and the grounded grid between the suppressor 
grid and collector plate (proportional to the yield at energy 

E primary)' 

More importantly, secondary electrons can be produced at the 
modulator grid or support ring. In fact, given the relatively 
large area of the modulator support ring compared to the area of 
the grid, most of the secondaries presumably would originate on 
the inner surface of the support structure. Here, of course, the 
underside of the grid cannot be illuminated by primary electrons 
as it can by reflected photons. The current so produced will be 
proportional to the yield at the energy E primary eV'morf- The 

total cuirent so produced is actually given by (cf. eqn. 13) 


= «J dv, v, J f d*v,/(v) Tb (-j- mv 2 - eV„J) (28) 

w m ""** ""** 

This current will flow with an annular cross section (cf. Figure 
94), intersect the suppressor with an energy of E vrvwy - eV^ 
(recall that the secondaries are "bora” with energies on the order 
of 1 eV), and continue to the collector plate with this 
intermediate population of "daughter" secondaries. The current 
reaching the collector plate can be modulated due to i) 
modulation of the primary flux reaching the modulator grid 
and/or support ring, ii) modulation of the yield at the modulator 
grid and/or support ring, iii) modulation of the number of 
electrons getting by the shield ring, and iv) modulation of the 
yield at the suppressor grid. In regard to mechanism iii), some 
electrons may be stopped by the shield ring. An assessment of 
possible modulation due to the angular distribution of electrons 
emitted from the ring support combined with the potential 
distribution above the suppressor grid is beyond the scope of 
this study, however. 

Figure 96 also shows possible trajectories of electrons 
produced by energetic primary electrons. Again, in this case it 
may be possible to produce a current on the collector plate 
which follows the yield curve for gold. This is apparently 
necessary to explain the "turn on" in the E2 plasma data 
observed in at least one of the "voids" observed by Voyager 2 
near Ganymede. 

The only source of readily available data on the relative 
amount of modulation in the secondary electron current for both 
normal and reversed suppressor configurations comes from early 
in the mission when both forward and reversed suppressor 
configurations were used in conjunction with DC return 
measurements. There are also a few places where both DC 
return and plasma data are available for the case of large 
secondary currents produced by energetic elections. These 
cases occurred in the Jovian magnetosphere during the Voyager 
1 pass. The "clean" data, i.e., those in which the fluxes are 
high enough to be over the DC threshold and positive ions are 
not flowing into the D cup, occurred outbound on day 65 of 
1979. 

First we consider spectra in the solar wind near the earth 
when the ambient electron population was relatively cold, i.e., 
no ambient electrons were energetic enough to produce a real 
signal in the E2 mode. The following table lists measured E2 
currents as a function of channel number using the DN to DC 
return conversion of Table 1. All currents are in units of fA 
(10 ~ 15 amperes), and TH indicates a threshold current for the 
measurement in question. These spectra are plotted in Figures 
97 through 100 along with neighboring El and L spectra (A and 
D cups). The format is the same as that used for Figures 55 
through 69, with the currents from the plasma measurements 
plotted as positive for both ions and electrons. These figures 
also show the appearance of the solar wind in the A cup during 
the DC return measurements. 
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Vojrtgv 2 1977 

Suppreaaor 


Normal 

Reveratd 

Time: 

243:1856 


238:2011 

243:1954 

242:2332 


l DC 


l PLS 

l DC 

tpLS 

1 

-1.00 X 10 s - 

TH 

8.475 x 10 2 

-137 x 10 5 

1.740 xlO 3 

2 

-1.00 x 10 5 - 

TH 

9.787 x 10 2 

-138 x 10 5 

23*5 xlO 3 

3 

-600 x 10 5 


1.353 x 10 2 - TH 

-8.00 x 10 5 

1 .353 xlO 2 - TH 

4 

-3.20 x 10 6 


1.353 x 10 2 - TH 

-4.48 xlO 6 

1.353 xlO 2 - TH 

S 

-5.44 x 100 


1.052 x 10 5 

-8.32 x ID 6 

1353 xlO 2 -™ 

6 

-3.92 x 10 6 


1.353 x 10 2 - TH 

-960 x 100 

1353 xlO 2 - TH 

7 

-6.40 x 100 


1.353 x 10 2 - TH 

-960 x 100 

1.353 xlO 2 - TH 

8 

-6.40 x 10 s 


1.353 x 10 2 - TH 

-960 x 100 

7.339 xlO 4 

9 

-6.40 x 1C0 


1.355 x 10 2 - TH 

-960 x 1C0 

2339 xlO 5 

10 

-6.40 x 100 


1.172 x 10 5 

-960 x 100 

437* xlO 5 

11 

-6.40 x 10 6 


2.494 x 10 5 

-*.96 x 100 

5.705 xlO 5 

12 

-5.92 x 100 


2.160 x 10 5 

-*.32 x 100 

5304 xlO 5 

13 

-5.92 x 100 


1.305 xlO 5 

-76* x 100 

437* xlO 5 

14 

-5.92 x 10 6 


4397 x 10 4 

-76* x 1C0 

2.986 xlO 5 

13 

-5.92 x 100 


1.353 x 10 2 - TH 

-7.04 x 100 

1 -939 xlO 5 

16 

-5.44 x 100 


1.353 x 10 2 « TH 

-704 x 100 

1353 x!0 1 -TH 


It is interesting to note that electrons produced at the 
modulator cannot reach the collector plate for the lowest 3 
channels due to the suppressor potential. The currents measured 
in those channels could be due to feedthrough from higher 
energy electrons, but these are not present if “high energy'’ 
electrons are simply those photoelectrons ejected from the 
modulator grid. A constant current in the reversed mode would 
simply be due to photoelectrons produced at the suppressor grid, 
but here that grid is not illuminated. At face value there seems 
to be no way of producing a negative current In these 
channels; this is one mystery we return to later. 

One final thought on this subject the angle between the 
normal and the direction to the sun is 73. 3°, as stated above. 
For sunlight not to be fully intercepted by the shield ring and 
the structure of the cup above it, an angle of less than 70.2° is 
required. Obviously, sunlight cannot be producing 
photoelectrons off of the suppressor grid, and these would not 
be modulated anyway. However, the solar wind is entering the 
cup at a large angle of incidence, and, as the modulator 
potential is increased, the protons are refracted more toward the 
center of the cup. This refraction effect may be just enough to 
illuminated the suppressor grid, produce secondary electrons via 
proton bombardment (more of this later), and produce a 
negative current from the secondaries. The protons themselves 
would have to miss the collector plate itself for the net current 
to be positive as the secondary yield is below unity for typical 
solar wind proton energies. In the corresponding positive ion 
mode, the DC return current in the side sensor is at positive 
threshold, while the measurements in the A cup suggest the 
solar wind current is -10* f A. With all of the efficiency factors 
involved, this mechanism appears to be ruled out, as negative 
current approaching this in magnitude is measured in the D cup. 

The solution to a complicated problem is always the obvious 
one! Recall from above that the measured suppressor potentials 
were -95.2 V in the normal mode and -85.9 V in the reversed 


mode. For channel 3, the modulator potential is oscillating 
between (nominal) potentials of 56.7 V and 92.3 V. Hence, we 
expect that a finite amount of secondary current can reach the 
collector in the reversed configuration. If the true modulator 
potential is off by as much as 3%, electrons would be just at the 
threshold for producing a current Currents in the first two 
channels in the reversed configuration are still an 
embarrassment 

A source of negative current in the third channel of E2 is also 
present during the Cruise Science Maneuvers (normal suppressor 
configuration). A spot check of the DC return data shows 
positive currents in the first three channels of E2 when the solar 
wind is entering the D cup; channels 4 through 16 have negative 
currents due to the larger currents produced by photoelectrons. 
However, while the currents in the first two channels are the 
same, the current in channel 3 is consistently lower. For 
example, for Voyager 1 at 1978-258/0139: 8.237 the angle to 
the sun is 43.2° (the minimum reached). The previous L mode 
has a DN of 153 in channels 7 (solar wind + suppressor 
inefficiency + photoelectrons from suppressor grid). The E2 
mode has DNs of 154 in channels 1 and 2, DN of 148 in 
channel 3, and DNs ranging between 0 and 34 (net negative 
current) in all the other channels. It is curious to note that the 
DNs in the previous El mode, all due to the solar wind plus 
photoelectrons, range from 165 to 174, i.e., much more positive 
current when the suppressor potential is lower Using the 
nominal values of Table 1, channels 1 and 2 of E2 show 
currents of 2.75 x 10 s fA, while channel 3 has 2.17 x 10 s fA, a 
change of 0.21. Given the relatively large suppressor grid 
transparency, this cannot be simply accounted for by secondary 
emission. Further study is required to pin down the source of 
negative current Another possibility is leakage due to the grid 
plane not being an equipotential (A. J. Lazarus, private 
communication); it is possible that this is the source of the 
current (from electrons produced on the modulator at 92.3 
volts). This mechanism cannot, of course, be the source of the 
negative current measured in the first two channels in the 
reversed suppressor configuration. 

As noted above, the normal suppressor DC currents are, at 
most, one third of the current measured in the reversed 
suppressor mode. The normal plasma measurements are, at 
most, -40% of the reversed suppressor plasma measurements, 
although this percentage is a fairly strong function of channel 
number. 

One can also consider the fractional modulation by by taking 
the ratio of the plasma to DC return currents, keeping in mind 
the possible uncertainties in the DC currents due to temperature 
dependence and the quantization of the steps (-10% as 
compared to less than 1% for the plasma measurements). These 
ratios are shown in the following table for those channels not 
registering only threshold current; for the latter N/A (not 
available) is indicated. 
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Voyager 2 1977 

Suppressor 

Normal 

Reversed 

Time: 

243:1856, 238:2011 

243:1954, 242:2332 

Channel 

- IplsHdc 

- Ipts / Idc 

r- 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 


N/A 

N/A 

N/A 

N/A 

0.01934 

N/A 

N/A 

N/A 

N/A 

0.01831 

0.03897 

0.03649 

0.02204 

0.00776 

N/A 

N/A 


0.01270 

0.0202 

N/A 

N/A 

N/A 

N/A 

N/A 

0.00764 

0.02332 

0.04456 

0.06367 

0.06615 

0.05570 

0.03888 

0.02574 

N/A 


We estimate that AS = 0.1 across one E2 channel near the 
peak plasma current, where A5 is the variation in secondary 
yield (cf. eqn. (10)). If the yield variation were the only source 
of modulation for 8 = 1, we would obtain 

/oc = /«c<5> < 29 > 

and 

Ins = /« AS [1 - (1 - ■^) 1 ] (30) 

which yields a modulated fraction I ms lloc x 10" 5 , a smaller 
value than is actually measured in the spectra considered above. 
We have assumed that the conditions producing the spectra are 
constant and that eqn. (10) adequately represents the yield. 
These are probably not bad assumptions. It is also possible that 
the DC return currents are larger than assumed above, i.e., that 
the fractional modulation is smaller than indicated in the above 
tables. Given that table 1 is for low temperatures and the 
thermal dependence shown in Figure 9, this is likely to be part 
of the problem. More detailed calculations of secondary 
electron trajectories within the cup and/or actual measurements 
with the flight spare instrument will be required to resolve the 
apparent discrepancy. 

As we have noted, a similar type of analysis is possible in die 
case that we know photoelectrons are not responsible by using 
data obtained during the Voyager 1 encounter with Jupiter. 
Currents and the modulation ratio from an example are listed in 
the following table (in the normal suppressor configuration). 

The fractional modulation is apparently comparable to that 
obtained when photoelectrons are present The variation with 
modulator potential is, of course, different in the two cases, and 
for both cases time aliasing, and the conversion of DC return 
DNs to current are possible sources of error. However, it 


appears that the resulting E2 spectra arc qualitatively consistent 


Voyijtr l 1979. D«y 63 

Chwnel Number 

DC Return Current (fA) 

Plum* Car mu (fA) 

Fraction*! 


1122:35.4)6 

1121:27.360 

Modulation 

1 

-1.09* 10 s 

1 -647 x 10* - TH 

N/A 

2 

-1.12 x 10 s 

1.647 x 10 2 - TH 

N/A 

3 

-120* 10 5 

1.647 x 10* - TH 

N/A 

4 

-1.28 k 10 5 

1.647 k 10* - TH 

N/A 

3 

-1.40 k 10 5 

1 ,647 k 10* -TH 

N/A 

6 

-1.40 k 10 5 

1.647 x 10* - TH 

N/A 

7 

-l .40 k 10 5 

1047 k 10* -TH 

N/A 

s 

-1.45x 10 5 

3.307 x 10* 

0.0024 

9 

-1.40 k 10 5 

8011 k 10* 

0.0062 

10 

-1.40 k 10 5 

1.708 x 10 3 

0.0122 

11 

-1.40 k 10 5 

2.726 x JO 3 

0.0195 

12 

-1.36 k 10 5 

3007 x 10 3 

0.0238 

n 

-1.21k 10 5 

3.264 x 10 3 

0.0255 

14 

-120* 10 5 

3.769 x 10 3 

0.0314 

15 

-103 x to 5 

5.210 x 10 3 

0.0496 

16 

-100 k io 5 - th 

5,401 x 10 3 

N/A 


We have also examined the El spectra at this time. The DC 
return spectrum at 65/1123:43.416 exhibits a broad maximum 
from channel 4 through channel 8; the plasma measurement at 
65/1122:15360 shows occasional currents above the 0 DN level 
which have no organized structure, i.e., they are consistent with 
noise. 

Given that there were no DC return measurements in the the 
Voyager 2 measurement sequence at Jupiter near the so-called 
voids, but there is also no indication of random noise in the 
lower channels, i.e., the lower channels appear to be "clamped" 
to threshold values, we con only conclude that those spectra are 
produced by a very hot primary population which produces 
significant numbers of secondary electrons within the instrument 
housing. 

The El Spectra 

We have not concentrated on the El spectra for because the 
physics of secondary emission is more complicated for 
impacting electrons at under 50 eV. In particular, it is difficult 
to differentiate between secondary electron, backscattered 
electrons, and reflected primaries at these low energies. It is 
clear from Figures 55 through 58, 97, and 99 that in the El 
spectra, the current at the collector increases for increasing 
potential for at least the first two, and sometimes, the first three 
channels (typically for the reversed suppressor configuration). 
This behavior covers the range from 10 to ~25 volts. However, 
from -25 to - 140 volts, the measured current decreases with 
increasing modulator potential, i.e., the collected current does 
not simply follow the expected yield curve as occurs in the E2 
mode measurements. This behavior is clear in the PLOMAN 2 
spectrum from 1979-62/1832:31.471 as well (angle to the sun 
43.7°). In the cruise science maneuvers, the signature is 
complicated by the solar wind at such small angles: the total 
current in the El spectra are positive in all channels but with 
variations due to the ambient electrons and photoelectrons also 
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p resent. In addition, the secondary yield is expected to be 
small; it is less than unity for less than 250 eV electrons 
bombarding tungsten [e.g., McDaniel , 1964]. 

The other fact of interest to note is that on the figures referred 
to immediately above, the £1 and £2 spectra do correspond. In 
the normal suppressor configuration, -the same current is 
measured in the fourth channel of £2 and the top four channels 
of £1. In the reversed configuration, the current measured in 
£1 is always somewhat larger in magnitude. This variation 
between the two configurations suggests part of this behavior 
reflects the inefficiency of the suppressor in the normal 
configuration. 

The decrease in measured current with increasing energy 
above 25 volts on the modulator must, therefore be related to 
the suppressor potential, by elimination. One possible way of 
decreasing the negative current would be by increasing the 
energy of the electrons liberated as secondaries from the 
collector plate as the primary energy increases. True 
secondaries have a common energy distribution independent of 
the energy of the incident primary. At 25 volts, a substantial 
fraction of secondaries from the collector can pass the 
suppressor to be absorbed at other grounded surfaces in the cup. 
The fraction lost is about 8% for the reversed suppressor and 
and 10% for the normal configuration (see Appendix B and 
Figure 15). 

The energy of the secondaries is not a function of primary 
energy, but the energies of backscattered electrons are 
([Sternglass, 1954] and references therein). Unfortunately, the 
literature on secondary and backscattered electrons at energies 
below 200 eV is limited and probably dated with respect to the 
measurement techniques employed. It is interesting to note that 
there are secondary maxima in the energy spec eras of secondary 
electrons which are characteristic of the target material, although 
this seems to be related to the prior heating of the target 
[Hacheriberg and Brauer , 1959]. Rudberg [1936] found such 
secondary maxima for a gold target at -3 and -24 eV. If such 
maxima are present in the secondaries emitted from the collector 
plate, then as the modulator potential (at which the 
photoelectrons are produced) exceeds 24 volts, electrons at this 
characteristic energy can be emitted; they would immediately be 
lost as they could not be stopped by the -8 volts on the 
suppressor grid. If a larger fraction of the secondaries are 
emitted with this energy as the primary energy increases, this 
could account for the decrease in negative current at higher 
energies in the £1 mode. However, this does not appear to be 
easily reconciled with similar measured currents at similar 
accelerating potentials in the £2 mode; with a suppressor 
potential of -95 V, all the the higher energy secondaries should 
be stopped. This change in suppressor efficiency should be 
reflected in an abrupt increase in collector current magnitude. 
Such a change is, of course, not present 

A comparison of the £1 and £2 spectra does suggest that 
whatever mechanism is responsible for electrons to leak past the 
suppressor to the collector plate in the E2 mode is also 
responsible for them leaking past the suppressor in the El 
mode. This process must become independent of the suppressor 
potential once the primary energy has exceeded -90 volts. 
Again, this could be due to the details of the secondary, 


reflected primary, and backscattered primary electrons at the 
energies where the distinction between these various populations 
becomes blurred. 

A final consideration for interpreting El spectra is a problem 
with the amplitude in the £1 measurements (see section 3.1 of 
report by R. J. Butler, unpublished manuscript). The effect is 
present at some level in the first two channels and is due to 
transients from switching the suppressor potential just prior to 
the El measurements. Calibration curves exist which can be 
used to remove this problem in the analysis of the electron data. 

Secondary Electrons Produced By Positive Ions 

A subject not touched upon here, but of potential importance 
in the analysis of some of the magnetospheric data obtained by 
the Voyagers, is the production of secondary electrons by 
positive ions. Production of such electrons at the modulator grid 
was of concern in the analysis of IMP-2 and OGO 1 and 3 data 
obtained by M.I.T. Faraday cups in the Earth’s magnetosphere 
[ Vasyliunas , 1968; Binsack, 1967]. Vasyliunas [1966] noted 
that due to the energy dependence of the secondary yield with 
respect to ion energy, the current produced will be out of phase 
with respect to the synchronous detector (this effect has also 
been noted by Sittler [1983]). As the yield is an increasing 
function of proton energy, more negative modulator potentials, 
should lead to larger yields of electron currents. Evidence that 
this was indeed occurring came in the form of electron data 
from OGO 1 in which the telemetered signal was clamped to 
the instrument zero level, below the noise level [Binsack, 1967]. 
These measurements occurred in the plasmapause where all of 
the cold, high density protons appeared within a given energy 
window, exceeding the nominal electron flux in that window by 
some two orders of magnitude. 

The behavior of the measurements, i.e., a "clamping" at the 
zero DN level, is exactly the phenomenon found in the E2 
measurements discussed above. In that case, the electrons have 
to simply be "hot", i.e., mostly in an energy range over which 
the secondary yield coefficient exceeds unity. 

The yield coefficient for 0 + incident on tungsten exceeds 
unity around a few keV [Large, 1963]. It is almost certainly 
the case that that the 0 DN values in the upper channels of the 
E2 spectra in the warm Io torus (1979, day 64 for Voyager 1) 
are produced by this effect, i.e., the positive ions are accelerated 
to -5 keV by the negative modulator potential, and hit the 
modulator grid yielding secondary electrons. This current 
would be out of phase with that due to the primary electron 
population, and, if of sufficiently large magnitude, can produce 
the observed effect [see also Sittler and Strobel , 1987], A more 
complex version of this effect, including yield curve modulation 
of the daughter secondaries, is almost certainly responsible for 
the bizarre E2 spectra recorded in the cold Io torus when large 
fluxes of positive ions could reach the (negative) modulator 
support ring but not the collector plate itself (due to the 
geometry of the flow). 

Comparison with Other Instrumentation 

The problem of secondary electron contamination in plasma 
experiment observations in space is not limited to modulated 
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grid Faraday cups. Difficulties associated with secondaries in 
RPAs due to large positive ion fluxes [Knudsen and Harris, 
1973] as well as in electrostatic analyzers due to large electron 
fluxes [ZwicU et al., 1987] have been reported. In both cases, 
the problem was only recognized when trying to explain 
unexpected observations. This is also true of the motivation for 
this report; however, our investigations indicate that the effects 
of secondaries can be present at other times as well and not be 
noticeable due to the "masking" effect of expected observations. 
Such behavior is, in principle, possible in RPAs and 
electrostatic analyzers as well, especially in exotic environments 
for which only limited data sets are available. Hence, one is left 
with the caveat that significant corrections for spurious 
secondary electrons may be required to characterize the 
local pi°«n» environment around a spacecraft, but the need 
to make such corrections may not be obvious in the data. 

Remaining Questions 
Quantitative Modeling 

The central question is whether a quantitative model of these 
secondary emission effects can be made which will allow us to 
simulate the data to a high degree of accuracy (-1% level). 
Such simulations are possible with the positive ion date in all 
four Faraday cups under conditions of negligible charging in the 
Jovian magnetosphere even when the ions are hot (see, e.g., 
[Barnett, 1986]. To accomplish this task the cascading or 
"feedthrough" of higher energy ions into lower energy 
measurement channels must be effectively taken into account. 
The feedthrough of electrons into the positive ion mode and of 
positive ions into the electron mode, both of which have been 
known of for some time [Vasyliunas, 1971], have been 
incorporated into the simulation code PMODEL (written by V. 
M. Vasyliunas, with documentation and corrections by J. W. 
Belcher). As noted above, PMODEL has been modified (with 
the efforts of G. S. Gordon, Jr.) to include the effect of 
suppressor inefficiency in simulating electron spectra. The 
ii ^Hitinniil current source provided by secondary electrons 
created at the modulator grid/support assembly and the 
suppressor grid have not as yet been included. 

The possibility of feedthrough from higher to lower energies 
was also recognized by Vasyliunas [1971] but was dismissed as 
a "higher order correction". The latter is certainly the case for 
measurements of the supersonic solar wind, the task for which 
the PLS experiment was primarily designed [Bridge et al., 
1977], but is a crucial effect to model in interpreting plasma 
data returned from the magnetospheres of the outer planets. 
Such a feedthrough effect is also present in electron data but is 
not modeled in the current analysis scheme (details of the 
analysis are given by Sittler [1983]). As more superthermal 
electrons are added this effect becomes more important; it 
accounts for -10% corrections to the relatively well-behaved 
electron data acquired during the Uranus encounter of Voyager 
2 (see Voyager Memorandum 136). A recent analysis of the 
electron date from the Io torus by Sittler and Strobe! [1987] 
showed that ion feedthrough provides a substantial contribution 
to the high energy electron channels. The calculation was made 
with the unmodified PMODEL code, i.e., suppressor inefficiency 


not taken into account, and was based upon the thermal ions; 
the feedthrough effect is actually larger as there are significant 
fluxes of non-Maxwellian ions in this region [ Bagenal , 1987], 

The other effects of secondaries with which we are concerned 
have been recognized in a qualitative way for some time as 
"extraneous effects" characteristic of the modulated-grid Faraday 
cup scheme [ Vasyliunas , 1971]. These effects include the 
modulated signal produced when sunlight illuminates the interior 
of the cup, which was not well understood at the time but 
thought to be related to yield of secondaries produced at the 
suppressor by photoelectrons produced at the modulator grid; 
this mechanism appears to be solidly supported by the Voyager 
data. They also include the out of phase signal resulting from 
secondaries produced by ion impact on the modulator grid (in 
the electron mode), and the enhancement of electron current due 
to secondaries produced at the modulator and suppressor grids 
in the electron mode. That the latter effect can also produce out 
of phase signals was not noted or, more likely, not considered, 
before as a very hot primary election population is required to 
produce this effect Such a situation has been encountered in 
the magnetospheres of Jupiter and Uranus (see below), and the 
study of these bizarre plasma conditions provided the impetus 
for this in depth study. 

The modeling problem is also complicated by having no a 
priori information as to the functional form of the distribution 
function to be used. The current state of affairs is illustrated in 
Figure 101. We show El and E2 spectra taken in the Jovian 
magnetosphere at the top. Note the 0 DN levels in channel 1 of 
El and channel 3 of E2. Both of these features, as well as 
potentially low values of current in channels 2 and 3 of El and 
4 of E2 are due to out-of-phase signals produced by internal 
secondary emission due to the electrons present in the non- 
Maxwellian tail. The second two plots (from the top) show a 
best simulation from the unmodified PMODEL code. Note that 
there can be no direct current in the first 3 channels of E2 due 
to the suppressor potential; the currents there are entirely due to 
feedthrough from higher energy electrons. Note also that the 
simulated and measured currents in the first two channels of E2 
show a rather good match. This simulation assumed a k 
distribution (see , e.g., [ Vasyliunas , 1971] and references 
therein), with k = 1.5. The next two rows of plots show, 
respectively, a simulation with k = 5.5 and k — ► i.e., a 

Maxwellian, both with only somewhat higher most probable 
speeds. Finally, the fourth row shows the simulated two 
Maxwellian fit from the NSSDC. Although, the El currents 
match fairly well, the match to the high energy tail is not good. 
This latter discrepancy has been alleviated to some extent by the 
use of four Maxwellians to fit the data (E. C. Sittler, Jr., private 
communication). Note that none of these simulations include 
effects due to secondaries or ion feedthrough. Note also that 
the effects in the first few channels of the El spectrum 
introduce difficulties in the determination of an electron 
"temperature" (or most probable speed). This problem 
systematically occurs when the apparent electron temperature 
exceeds a few tens of volts (E. C. Sittler, Jr., private 
communication). 
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Other Qualitative Problems 

There is a great deal of data, especially from the Voyager 
encounters with Jupiter, which has not been analyzed because 
we have not understood how to do the analysis. The example 
of the "voids" near Ganymede during the Voyager 2 encounter 
is one, but not the only, example. The positive ion data 
acquired during the times in which 0 DN levels are recorded in 
the E2 data have been mystifying as well. The Voyager 2 
encounter with Uranus may have provided the Rosetta stone 
with which these data may be understood, at least qualitatively. 

Spacecraft Charging at Uranus . Figure 102 shows positive 
ion and electron data (long and short integration times) from the 
Uranus encounter during that time that a large (- -400 V) 
potential on the spacecraft is inferred [McNutt et al 1987]. It 
is important to note that the low energy tail is not "real", i.e., 
the ions are all accelerated to a potential consistent with the 
peak location in the ion spectra; the "tail" is purely a feed 
through effect due to these ions down to lower channels. The 
plot format is described in Voyager Memorandum 122 and is a 
modification of the "1800" plots produced during the other 
encounters (refer to Voyager Memorandum 79) to incorporate 
all of the GS-5 mode data. The charging signature is readily 
apparent in the D cup in all of the L and M mode spectra. 
Similar signatures can be seen in the A, B, and C cups. It 
should also be noted that the E2 spectra exhibit the "turn on" 
signature characteristic of secondary electron production. In 
other words, if the ambient electrons are hot enough to 
produce secondary electrons which contaminate the E2 
spectra, then the ambient electron flux Is large enough to 
charge the spacecraft to negative potentials of several 
hundred volts. The transition to clear signatures in E2 occur as 
the ion spectra show an average acceleration through -250 eV. 
The channel corresponding to the turn on is about one channel 
lower than that containing the the peak in the L mode. This is 
the qualitative behavior expected: the spacecraft charges to a 
potential corresponding to the peak location in the electron flux, 
and, if the latter is below the maximum in the secondary yield 
curve, the secondaries produced will provide a large out of 
phase signal in lower channels. Further quantitative analysis of 
these features may be possible; however, such an effort is 
hampered by the relatively small signal-to-noise ratio and 
additional noise produced by the LECP stepping motor. 

Inferred Charging at Jupiter, Figures 103 and 104 show -45 
minutes of data from the PLS experiment on Voyager 2 as it 
was outbound from Jupiter. The spacecraft is -21 R / from the 
center of the planet, the D cup should show some response to 
corotating flow of plasma ions, and the main sensor cups none 
at all. The obvious question is what is the origin of the currents 
in the positive ion measurements? The E2 spectra suggest an 
answer, in light of the Uranus charging data. In all of the E2 
spectra, channels 2 through 7 are clamped to 0 DN; in addition, 
the turn on occurs slowly with channel number suggesting the 
gradual reduction of an out-of-phase signal. Such electron 
spectra can only be produced if the ambient electrons are hot. 
The Uranus spectra suggest that electrons which are that hot 
provide fluxes sufficient to charge the spacecraft negative. 
Given the similarity of the ion spectra with those at Uranus, it is 


probable that the spacecraft is charged to a large negative 
potential, and most of the ion currents are the feedthrough 
signatures of the positive ions (of several different plasma 
species) which are being accelerated into the four cups of the 
PLS instrument The occurrence of a peak in the ion spectra in 
channels 14, 15, or 16 implies that the potential is at least -2500 
volts. There is roughly a 24 hour period during which these 
charging levels were maintained. If this interpretation is 
correct, it may be possible to extract positive ion densities 
which can be checked against the cut off frequency observed in 
the Plasma Wave (PWS) data [Scarf et al. t 1981]. 

Figures 105 and 106 show a similar set of spectra which 
include the last of the plasma depletions observed during the 
Voyager 2 inbound pass [McNutt et al ., 1987]. The E2 spectra 
show the same qualitative behavior but are clamped to 0 DN up 
to a larger channel number. Most of the positive ion spectra 
peak in channel 16 or higher (with the two obvious exceptions). 
Both of (he features are consistent with a hotter ambient 
electron population and, as a consequence, more negative 
spacecraft potential Such behavior is observed in all of the 
depletions. In all qualitative respects, the data obtained within 
the depletions is consistent with that obtained for more extended 
periods outbound illustrated in Figures 103 and 104. 

Similar data were obtained by Voyager 1 outbound from 
Jupiter. Figure 107 shows data obtained -22 R j from Jupiter by 
the PLS instrument on Voyager 1. At the top of the figure and 
just past the center are E2 spectra similar to those of Figure 
104. The corresponding L mode spectra also show the 
characteristic rise with channel number which would be 
associated with ions accelerated into the various Faraday cups, 
t the center of the figure are two sets of DC return spectra 
acquired -1120. One of these was used in the analysis of the 
fractional modulation due to secondaries given in tabular form 
above. It is readily apparent from the E2 spectra in Figure 107 
that the turn on in the plasma spectrum is associated with the 
location of the peak in the DC return E2 spectra. 

Puzzling Features in the Ion Spectra . The lower part of 
Figure 107 shows a feature which is not so readily understood. 
The E2 and El spectra indicate relatively cooler electron 
temperatures, and the lower channels in the ion spectra have 
apparently been filled in with direct signals. However, Channel 
16 in the main sensor cups remains at a relatively low level. 
Figure 108 shows similar behavior a few hours earlier on day 
65, along with some bizarre jumps in the data in the D cup (the 
glitches in the M mode spectra are from LECP 6 second 
stepping). Figures 109 and 110 show similar behavior from the 
inbound pass of Voyager 1 on day 63 of 1979. In this latter 
case, azimuthally flowing ions are apparent in the D cup. 

Similar E2 and L spectra were also acquired for a limited 
region during the Voyager 1 outbound pass at Saturn. The 
signatures are not as clear due to the signal to noise ratio. 

The ion spectra in Figure 110 are particularly puzzling in that 
there is a drop to the 0 DN level in high channels of the C cup 
for both the L and M modes, without corresponding behavior in 
the A or B cups. The 16th channel of E2 also seems to be 
adversely affected. These spectra are not easily understood. If 
the signals in the upper channels of the A, B, and C cups in the 
L mode were due entirely to the ambient plasma population, the 
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ions producing the signature would be hot, and the signature 
should be the same in all of the cups. Large anisotropies can 
also be ruled out if the anisotropy is gyrotropic. Similarly, if 
the measured currents were due to energetic electrons, a similar 
signature should be present in all the cups of the main sensor, 
which it clearly is not. In addition, the affected cup is not 
always C, i.e., the cup showing 0 DN in channel 16 changes, so 
an electronics failure can be ruled out. 

The 0 DN levels certainly suggest an out-of-phase signal 
which must somehow be associated with the direction of 
positive ion flow. We come to this latter conclusion as there is 
no other phenomenon to distinguish between cups at the highest 
energies. Unfortunately, there are no DC return measurements 
for which the downturn is seen in the A cup plasma 
measurements and the fluxes are sufficiently above the 
photoelectron current (from suppressor leakage) for the change 
in net phase to be directly measurable. 

An out of phase signal is not readily produced in the positive 
ion measurements. As the modulator has a positive potential, 
secondary electrons can only be produced at the negative 
suppressor grid. For such a current to be modulated, the rate 
must be varied via the refraction modulation of the particles 
reaching the suppressor grid [Vasyliunas, 1971]. To produce an 
out-of-phase current, the net (positive) current must increase 
with incre asing (positive) modulator potential. Either, more 
positive ions must reach the collector plate or less electrons. 
Neither of these options has an obvious mechanism. The out- 
of-phase signal does seem to occur in the cups furthest away 
from the ion flow. As secondary electrons can escape only 
from grids charged to negative potentials, and the energies are 
two low to produce ions from either electron or ion impact, it 
would appear that the effect must be related to the collection of 
fewer secondary electrons with increasing modulator potential. 
Further study will require a full accounting of all feedthrough 
and secondary production corrections. 

Figures 111 and 112 show a sequence from Voyager 1 
inbound as the spacecraft neared the warm Io torus. Toward the 
end of the sequence, there is a component of convective flow 
into the B, C, and D cups. Up until -0135, there is a relatively 
high energy background signal in all cups in the positive ion 
modes. At this time the superthennal electrons disappear and 
both the electron and ion spectra change qualitatively. The 
signals in the high energy ion channels drop, indicating that 
these signals were probably feedthrough from the hot electrons 
(those with energies greater than 95 volts which could reach the 
collector in spite of the suppressor potential). After this time 
much of the electron currents may be due to ion feedthrough in 
the upper channels; the 0 DN level in channel 16 is almost 
certainly due to secondary electrons produced by the ions at the 
modulator grid. This effect must also have been present in 
previously but masked by the hot ambient electrons. Such a 
combination of effects would explain the somewhat depressed 
electron current in channel 16 of the E2 spectra in Figures 103 
through 106 and 109 through 111. In some of these figures, the 
ions producing the secondaries must have been accelerated by 
the spacecraft potential in the first place where the latter was 
caused by the hot electron flux; such synergistic effects 
apparently result in the complexities found in attempting to 


explain many of the ion and election spectra. 

Other strange features are found in spectra obtained with the 
Voyager 1 instrument within the Io toms itself. Figures 113 
and 114 shown data from -9 R / and -6 R j from Jupiter, 
respectively. In Figure 113, out-of-phase currents contribute to 
the low energy channels in the A and B cups in the L mode but 
not the M mode. In Figure 114, intermediate channels of the L 
mode spectra from the B and C cup exhibit large subtracting 
currents never seen again in the mission. This behavior could 
be related to the saturation of the amplifiers in the L mode, as 
there is no analogous affect in the M mode spectra and the L 
trade measurements are highest in the B and C cups of the M 
mode. More thought should be given to this possibility. 

The Oscillator Effect . Figures 115 and 116 show the best 
example of the most bizarre effect to appear in the data. For 
this particular period, the E2 spectra are relatively unaffected, 
although there are counterexamples. The El, L, and, especially, 

M mode spectra exhibit oscillations with respect to channel 
number. In the M mode, the oscillation period is clearly 
resolved, with the telemetered currents reaching down to the 0 
DN level between peaks. This behavior is presumably to be 
interpreted as ind icting an out-of-phase signal in these regions; 
if DC return measurements were available (and none seem to 
be) they would apparently show a current varying sinusoidally 
with respect to the channel number. It is worth noting that the 
widths of the zero DN regions are less than half the widths of 
the nonzero DN regions. In other words, a the sinusoidal 
variation appears to be superimposed on top of an in-phase 
modulated current, shifting the oscillation point upward from the 
zero level. That this is likely the case is more evident in the 
bottom of Figure 116. The oscillation frequency appears to 
reach a minimum -192/0815, and a similar minimum occurred 
-192/0515. Plate 1 of Scarf et al . [1981] suggests that these 
times may be associated with rapidly changing plasma densities. 
During these times, the plasma frequency was rapidly swept 
through 2.4 kHz, which is the frequency of the square wave 
power supplied to the PLS experiment; the plasma frequency 
apparently remained above 1 kHz throughout this region. The 
square wave applied to the modulator grids has a frequency of 
400 Hz, so it is not obvious why an effect is seen. Presumably, 
the oscillations in the measured current are a beat phenomenon 
occurring at the difference between the local plasma frequency 
and some reference frequency produced by the PLS experiment. 
The 400 Hz square wave on the modulator grid will also 
produce odd harmonics, so the relevant frequency provided by 
the instrument may be 1.2 or 2.0 kHz (3rd or 5th harmonic, J. 
H. B insack, private communication). Only narrow band data 
from the PWS experiment is available at this time, so the error 
ban on the density determined from the cutoff of the continuum 
radiation are rather large. Further, more detailed comparisons 
are required to determine the exact limits which can be put on 
the density. 

A better determination of the relevant spacecraft frequency 
may be obtainable from Voyager 2 data in the solar wind just 
prior to its encounter with Saturn. Here the solar wind is 
clearly visible in the M mode spectra, as are the oscillations. A 
fit to the solar wind proton peak at 1981-236/1300:58 (J. D. 
Richardson, private communication) yields a proton density of 
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0.0475 cm' 3 . The oscillation envelope occupies -10 M mode 
channels, or -2.4 seconds. The corresponding electron plasma 
(from charge neutrality, neglecting the a particles) and 
difference frequencies would then be 1960 Hz and 0.42 Hz. 
This is suspiciously close to the 5th harmonic. The 3rd 
harmonic corresponds to an electron density of 0.018 cm -3 . 
Such low density regions have been sampled, and there is, 
presumably, more power in the third harmonic than the 5th. It 
is possible that such lower densities do not provide enough 
electrons to yield a signal above the noise level. It is also 
possible that the spectra of Figures 116 and 117 correspond to 
such lower densities. Again, careful comparisons are needed to 
further pin down this effect We note that a plasma frequency 
of 400 Hz corresponds to a density of 0.00198 cm' 3 . Rather 
large speeds would be required to produce a measurable flux. 

For the solar wind spectrum considered above, the amplitude 
corresponds to a flux of -2 x 10 s cm' 2 s* 1 . For the density 
quoted above, the speed is 4.2 x 10 6 cm s’ 1 . This does not 
seem to relate to much of anything but is based on the 
assumption that the oscillating charges are spread out uniformly 
on the collector. If the area of the collector impacted were 
smaller, the corresponding particle speeds would be larger. 

Obviously, this phenomenon is not understood. Other 
complications include the exact mechanism of setting up current 
oscillations at the collector plate and to what extent, if any, 
Doppler shifts in the plasma frequency affect the exact 
resonance conditions. The reason for selecting a 400 Hz 
frequency for the modulator was a concern about resonance 
effects (A. J. Lazarus and J. H. Binsack, private 
communication). The alternating changes in the modulator 
potential are certainly suspect for establishing charge density 
oscillations, but no detailed studies of mechanisms were ever 
conducted. Another possible clue is the fact that the FWS 
experiment M s ees" the PLS experiment when the modulator is 
on; presumably this is due to waves excited in the ambient 
plasma (H. S. Bridge, private communication). 

Implications 

There several implications which follow from this study, one 
being major. The major result is that the energetic electron 
spectra, i.e., those obtained in the E2 mode cannot be taken at 
face value. We have also noted that corrections to the £1 
spectra can be important if the electron temperature is high 
enough. Even in those cases where there is no obvious out-of- 
phase current, the electron current below channel 8 will be 
decreased and above channel 8 increased over the fluxes present 
in the parent population outside the instrument housing. 
Obviously, the extent of the corrections are a function of the 
electron distribution function which cannot be simply read off of 
the measured currents due to the electron cascade-type of 
feedthrough effect Further complications arise due to 
feedthrough from the ions, which can be affected by 
feedthrough from the electrons, and spacecraft charging. These 
effects, and the resulting corrections to the analysis, are 
certainly minimal in the solar wind or cold Io torus. In other 
locales they may also be unimportant to some level, e.g., 1% or 
10% say, but this must be judged on a case by case basis. 


These complications may (or may not) have a major effect cm 
the electron impact ionization rates for multiply charged ions in 
the warm Io torus as calculated by Sutler and Strobel [1987]. 
Their rates are probably within a factor of 3 or so of the true 
rates, but a more detailed study is necessary to confirm or reject 
such a number. Given that the ionization cross sections 
typically peak around 100 eV and secondaries apparently cause 
an out-of-phase signal up to -750 eV, the published rates are 
almost certainly low. Their other conclusions based upon 
details of variations in the hot electron component, rather than 
those based upon gross morphological characteristics, are also 
called into question by these results and should be reexamined. 

Consequences for positive ion results tend to be less dramatic, 
as a ready source of modulated secondaries within the 
instrument is no longer available. The major source of 
"contamination'' is expected to be simple electron feedthrough, 
which may explain currents in the main sensor in the upper 
channels, such as depicted in the top of Figure 111. The 
feedthrough contribution should be relatively flat (based on 
some simulations), and certain features remain unexplained, 
notably the downturn in the highest energy channels, present in 
the C cup in Figure 110. As a matter of course, we have 
refrained from publishing these data or any analysis based upon 
them because we could not duplicate them from our knowledge 
of the instrument response. This "policy" has, perhaps been 
overdone. A search far a qualitative mechanism for producing 
an out-of-phase signal in the ion measurements is required and, 
can perhaps lead to new physical insights concerning the 
characteristics of the ambient plasma. 

A final major consequence is that many of the spectra from 
the various planetary encounters can be understood in terms of a 
very hot electron environment and a negative charging of the 
Voyager spacecraft The work of Khurana et al. [1987] not 
withstanding, the spacecraft potential never seems to have 
exceeded significantly a few kV. Some of the spectra shown in 
Figure 105 appear at first sight to be difficult to reconcile with 
this solution, as a turnover in channel 16 alone must be viewed 
with suspicion. On the other hand, for the apparent feedthrough 
signals to be the amplitude they are with the corresponding 
spectral shape, the ion energy cannot appreciably exceed 5 kV. 
Some preliminary work suggests 7 kV is a hard upper limit 
unless appreciable plasma injection (factors of 10) is occurring 
simultaneously with the charging process. The latter can be 
ruled out as it would give rise to a corresponding discharge 
current to the spacecraft, lowering the potential, and 
contradicting the observations. The charging levels found in the 
Voyager data should prove useful in the interpretation of plasma 
data returned from the Jovian magnetosphere from the Galileo 
orbiter, as such charging levels had not previously been inferred 
from Voyager and Pioneer data sets (F. L. Scarf seems to be 
vindicated on this point!). 

Finally, we note in passing that the DC return measurements 
have been crucial in sorting out these various phenomena. At 
the same time, they are difficult to work with due to the 
(presumed) thermal dependence of the DN to current 
conversion. Examination of data during the thermal cycling of 
the Voyager 1 instrument seems to preclude as strong a 
dependence as inferred from the pre-flight calibrations, however 
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(see Voyager Memorandum #158). The IX return 
measurements are typically difficult to interpret due to the the 
presence of the suppressor grid and the combined effects of ion, 
electron, and secondary electron currents. The PLS experiment 
is fundamentally different from an RPA (retarding potential 
analyzer) due to the grounded aperture grid, i.e. f the instrument 
potentials do not, except to low order, intrude into the local 
plasma (this may be a bit misleading: using the formulas of 
Morse and Feshbach [1953], I estimate about a 4% "leakage” in 
the ion modes, or about 100 V for a modulator potential of 
2500 V). As a result, the floating potential of the spacecraft 
cannot be determined, and, at zero modulator potential, there 
should be a net zero current to the collector plate. That this 
does not happen in the Jovian and Saturnian magnetospheres 
can usually be traced to i) directional effects due to supersonic 
ion flow, ii) photoelectron production and suppressor leakage, 
iii) changing electron temperature with respect to the suppressor 
potentials used for El, E2, and ion measurements, and iv) 
secondary electron production at the skin of the spacecraft itself. 
Effect ii) accounts for the difference in the A and D cup 
currents in Figures 55 through 59. Effect i) accounts for the 
differences between the currents in the A and D cups in Figures 
61, 64, 65, and 66. Effects iii) and iv) apparently account for 
the zero in the current (with respect to modulator potential) in 
Figures 67, 68, and 71. It may be possible to understand the 
some aspects of the secondary emission characteristics of the 
spacecraft from these data, but it will not be an easy 
undertaking. Note that due to effect iii) the measured DC return 
current is almost always positive ; in the E2, L, and M modes, 
the suppressor potential is -95.2 volts, and current carrying 
electrons with energies less than this cannot reach the collector 
plate to cancel the positive current earned by the ions. This 
also explains why the D cup, L mode and E2 traces arc usually 
continuous as the average modulator potential changes sign. 
The spectra in Figure 61 arc the exception; in this case most of 
the electron current is carried by hot electrons. 

Summary and Conclusions 

We have completed an extensive survey of PLS data from the 
two Voyager spacecraft in an attempt to understand certain 
anomalous spectra obtained within the magnetospheres of the 
outer planets. The study of the plasma voids on the Voyager 2 
pass and the ensuing argument over charging levels with 
investigators studying LECP data motivated this study. In 
attempting to explain these spectra we have developed a much 
improved understanding of how modulated -grid Faraday cups 
actually respond to the often bizarre plasma environments 
encountered in deep space. This has included an understanding 
of the generation of photo- and secondary electrons within the 
cup housing, and how these populations affect the currents 
actually measured by the instrument and telemetered to Earth. 
Some of this understanding is quantitative while much is not 
Due to the various synergistic effects, such quantitative study 
will require the use of the flight spare in laboratory tests 
coupled with further numerical modeling. We have also found 
regions in these magnetospheres where the spacecraft was 
apparently charged negatively to kilovolt potentials. This may 


well open up possibilities for further analysis of PLS data. 
Finally, we have made the first systematic study of the DC 
return spectra from then PLS instrument, used them in 
interpreting some of the strange plasma data, and identified the 
difficulties and potential uses of this data set. 
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4.2108401 8.3820401 1.1109402 4.5709402 0. 9458402 8.4420402 1.0508403 1 1708403 

1.0009403 0.5408402 5.7860402 2.7180402 1.0508402 2.4900401 2.67Be400 1.0828-01 

TYPE 1 TO CONTINUE 
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I qr . tin. 1 i 

sA T iikYTri M 




isst :srs skis sss sirs ass ascs sscs 

9 VOYAGER DATA - CURRENTS (FA) 3- 4 

?:sc£ :.sss Mrs; irrss is 5:isr,; ises ;srs 

TYPE 1 TO CONTINUE 


3PL0T SIMULATED DATA? 1-YES, 0-NO 

She SCALE STARTS AT 1.00a+01 FEMTOAMPS AND GOES UP 3 DECADES 
0CHANGE CURRENT SCALE FROM CURRENT VALUE? 1-YES, 0-NO 


Plot option [cefghnqu]: g 


npstar sp 


1n3 gag [16 
1n3 gsg[17 
1n3 gsg[l8 
1n3 gsg[19 



SIMULATION ; 1S0TBPIC 
0.000 > KM/S B 


V = <0.000 , 0.000 , 0.000 > KM/S 
R * 0 Z * -1 

HPflfl'-11861?120 HPER - 11861.120 KM/S 


0.000 > COMMA 



Figure 46 






L T ii/ATHd:aiTVl 




0 . 0800408 0 . 000*400 0 . 000*400 0 . 000*400 0 . 000*400 0 , 000*400 0 . 000*400 0 . 000*400 

0 . 000*400 0 . 000*400 0 . 000*400 0 . 000*400 0 . 000*400 0 . 000*400 0 . 000*400 0 . 000*400 

3 VOYAGER DATA - CURRENTS (FA) J- 4 

3.590*401 5.460*401 9.5490401 4.0740+02 B.059©402 7.246b402 8.9760402 9.9858402 

9.1930402 7.3880402 5.077e402 2.4240402 9.4978401 2. 2880401 2.4730400 1.735e-01 

TYPE 1 A CONTINUE 


',Iih I ] tun I -,t 




0PLOT SIMULATED DATA? 1-YES, 0-NO 
1 

THE SCALE STARTS AT 1 . 00*401 FEMTOAMPS AND GOES UP 3 DECADES 
BCHANGE CURRENT SCALE FROM CURRENT VALUE7 1-YES, 0-NO 


Plot option [cefghnqw] ; g 


1n3 gsg[9] 


umpstar sp 



SIMULATION : COSINE 
Y = <0.000 , 0.000 , 0.000 > KH/S 
fl * 0 2 * -1 
N<CM-3> *.008 

WPflR -11861.120 HPER - 11861.120 KM/S 


B = <1.000 , 0.000 , 0.000 > CAMMP 


Figure 47 
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0 000*400 0 . 000*400 0.000*400 0 .000*400 0.000*400 0000*+00 ® ????**?? 

©I 000*400 0 . 000*400 0.000*400 0.000*400 0.000*400 0.000*400 0.000*400 0.000*400 


VOYAGER DATA 


CURRENTS (FA) 


4 401*402 7 . lfl7*+02 1.457*403 1.349*404 1.108*+04 8.047*403 7.403*403 0 252*403 
5 1321*403 41452*403 3.824*403 3.175*403 2.584*403 1.810*403 1.954*403 1.850*402 

TYPE 1 TO CONTINUE 



Figure 48 



firtSt /. a 


0.0009408 0.0009400 0.0009400 0.0009400 0.0009400 0.0009400 0.0009400 0.0009400 

0.0009400 0.0009400 0.0009400 0.0009400 0.0009400 0.0009400 0.0009400 0.0009400 

0 VOYAGER DATA - CURRENTS (FA) J« 4 

3.5039402 5.0739402 1.1739403 1.1009404 9.0909403 7.0309403 5.5449403 4.4929403 

3.7609^3 3.2179403 2.6329403 2.4299403 2. 8439403 1.5559403 6.8279402 1.5809402 

TYPE 1 TO CONTINUE 


0PLOT SIMULATED DATA? 1-YES, 0-NO 
1 

THE SCALE STARTS AT 1.009401 FEMTQAMPS AND GOES UP 4 DECADES 
9CHANGE CURRENT SCALE FROM CURRENT VALUE? 1-YES, 0-NO 


Plot option [c9fghnqu] : g 


FOR El MO 
IF (W 
CALCULATE 


999 FORMA 

|i 

15 RETUR 


END 

I 

0 

r 

BIMAXWELL 

r 




SIMULATION ; 

(90troplC 0099 

V * <0.000 , 0.000 , 

0.000 > KM/S 

B * <1.000 , 0.000 „ 0.000 > COMMA 

ft ■ 0 Z * -1 



N<CN-3> *.445 



MPflR -1861.019 NPEfl 

- 1061. 019 KM/S 



Figure 49 
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0/ WG«Yj! >; 

0f poor 


Quality 


B 000. 400 0000.400 0 000.400 0.006.400 0.0000+00 0.000-00 0- 000—00 0 000—00 

0.0000+00 0.0000400 0.000—00 0 . 000—00 0 000—00 0.000-00 0.000—00 0.000—00 

5 VOYAGER DATA - CURRENTS (FA) J- 4 

9 778040 1 1 5240402 2.0300402 1.8300403 1.9800403 2.0090403 1.9390403 1.7700403 

?:S42-03 l' 245— 03 oiSlSoK 7.4550402 5.414.402 3,554.402 1,716.402 2,781.401 

TYPE 1 TO CONTINUE 


r .h»* lltn.il -jt 


3PL0T SIMULATED DATA? 1-YES, 0-NO 

She scale starts at i.eea+ei femtoamps ano goes up 4 decades 

BCHANGE CURRENT SCALE FROM CURRENT VALUE? 1-YES, 0-NO 


Plot option [cafghnqui]: g 



SIMULATION iHaat isotropy* 

V * <0.000 . 0.000 , 0.000 > KM/S B = <1.000 

A * 0 l « -1 

HPPR -106KO19 WPER - L86L.019 KM/S 


0.000 > COMMA 



Figure 50 
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VOYAGER DC RETURN DATA 






V(j i HCE*n iDC'ntTUfin DAf 


^£X0 Z.8/0l/f>0 H OJta *»y aocd S J °i Ja > ua 3 'iTH '«T ’linNDW n -y 




VOYAGER DC RETURN DATA 


LZ 10 ^.8/OI/fiO H^esay aoedg joj ja>ua 3 -J.-PW -yp ‘imNDW *1 a 



Figure 
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ur 4 linN3W '1 *U 
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— . _J 
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VOYRGER DC RETURN DRTR 


oouu co/ui/nu 


jo*uaj ‘1 ' 1 'W 


’ur “iinNOw *i -y 





VOYRGER DC RETURN DRTR 


jd ^u/ui/nu 


1-juuwt.ay aoeog JOj JS^uaQ *j.*I*W 


•yr ‘iinNow n -y 








VOYAGER DC RETURN DATA 


X 1 w 


H! 11IINJW *1 -y 









VOYRCER DC RETURN DflTR 






VOYRCER DC RETURN ORTR 


Obd £hOO Z-8/OI/fiO 


i-jo-ieasay aoedg joj je>ueQ 


•«r *linN3W "1 'd 
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X SCALE 



Please hit RETURN to continue. (NOTE: Scree 







wish to continue drawing electric field lines 
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real(w) contours 0.100 
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Finure 


0 1255 — 0.4745 0.4520 
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contours O.lOOe— 01 
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real(w) contours 0.100 




0. 4530 . 
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Figure 83 
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Oreal(W) fgmtours 0.10fle+04 
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Oimag(W) .Contours 0.10Pe+04 



Oreal(Z) Contours 0.50Q£-01 




Oimag(Z) .Contours 0.104) 
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APPENDIX A 



Chg 7, 26 Oct 77 


¥ 


t 16-iQ J 

MJS 77 TELEMETRY CALIBRATION 


ft/3 aiihv 

PLS 


Calibrated or 

C, Odd 


ULASuRtMINT TlTlt 

PLS SEN 


SUBASSEMBLY NAME 

Plasma Sub system 


SINMNG BJ TNANSOUCU 
OIVICI _ 

□ OTHCB 


DATE 

10/21/77 


TTS NO 

E-622 


SUBASSY «tr NO 

SUBASSY SE* NO 

CAl (ORATION Tf MP 

2032 

002 

N/A 


TRANSDUCER type 

TRANSDUCER pah 1 NO 

IRANSOUtf* SC A NO 


ST 11785-0004 

8676 








OATC 


1 BAN50 EXCITATION 


VyHLI Ulxnl 1 L/ix • lu*u 
MEASURAND 

(units Leg. 'C 

OUTPUT 

(UNITS DN * 

ML ASURAND 

(units Deg. C 

OUTPUT 

'UNITS CN ) 

-157. 59 

83 

68. 299 

224 

-137. 50 

102 

89. 00 

232 

-116. 29 

120 

108. 19 

239 

- 96. 699 

135 

128. 50 

246 

- 75. 199 

150 

1 50. 19 

253 

- 55. 00 

163 



- 34. 70 

175 



- 14. 499 

186 



+ 5. 100 

196 



26. 199 

206 



46. 500 

215 




z 

Q 


co 


=) 

O 


260 
240 
2 20 
200 
180 
160 
140 
120 

100 

80 


.1 




TT~ 
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--l- 
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—U . 
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.. i. .. 
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♦ i 
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— r~ • 

i ■ ■ 


... y . 
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i — r 






— * — 1 — » 




-150 -100 

MEASURAND 


-50 


50 


(UNITS: 1( ^ c g- C 


150 
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E-622 


MJS 77 NO 
-2 

S/C-31 
VGR- 1 

JPL 3983 R 3/76 
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RUN 0697 0022PLS 


FRArt 15 





0 

2 

4 

6 

8 

— 



0 

-202.77 

-201.73 

-200.68 

-199.64 

-198.58 


MJS77 SERIAL NO MJS-2 

n_r<;77 faction pls sensor temp 


10 

-197.53 

-196.47 

-195.40 

-194.33 

-193 25 


FOS TREESW1TCH 10 .. 70 


20 

-192.16 

-191.06 

-189.96 

-188.84 

-187.72 


FDS RANCE 25-000 OHMS 

FOS SERIAL NO 2 


30 

-186.59 

-185.44 

-184.28 

-183.11 

-181.93 


SUBASSY REF NO. .... . 2032 


40 

-180.73 

-179.52 

-170.30 

-177.06 

-175.80 


SU0ASSY SERIAL NO... 002 






-170.60 

-169.25 



50 

-174.52 

-173.23 

-171.92 


MEAS CAL 10 RANGE ... “70 TO 100 DEC C 
ME AS CAL IB TEff* 


60 

-i67.ee 

-166.50 

-165.09 

-163.66 

-162.21 


XDUCER SERIAL NO 8676 


70 

-160.73 

-159.24 

-157.71 

-156.17 

-154.60 


XDUCER IMPEDANCE ... 

XDUCER CAL IB DATE .. 9 OCT 74 


80 

-153.00 

-151.37 

-149.72 

-148.04 

-146.34 


DATA PREPARED 0Y . . . ROSEMONT ENC CO 

D 

A 

T 

90 

-144.60 

-142.84 

-141.04 

-139 21 

-137.36 


APPROVED BY . . . 

APPROVED BY ... R. ROTTER 

100 

-135.47 

-133.55 

-131.59 

-129.60 

-127.58 


RUN DATE 0697 0022PLS 

A 

no 

-125.52 

-123.43 

-121.30 

-119.13 

-116.92 


STANDARD DEVIATION . .56643248-01 

N 

120 

-114.68 

-112.40 

-110.08 

-107.72 

-105.32 


COEFFICIENT AO -.20276987*03 

COEFFICIENT A1 52035174*00 

U 

M 

B 

130 

-102.08 

-100.39 

-97.87 

-95.30 

-92.69 


COEFFICIENT A2 25540733-03 

COEFFICIENT A3 12712726-04 

140 

-90.03 

-87.33 

-84.58 

-81.79 

-70.95 


MEASURAND OUTPUT F (X) 

E 

R 

150 

-76.06 

-73.13 

-70.15 

-67.12 

-64.04 


DN PRITC DELTA 

160 

-60.90 

-57.72 

-54.49 

-51.20 

-47.86 

— 

-70.00 154.15 -69.9296 -.070376 

-60.00 160.56 -60.0212 .021169 


170 

-44.47 

-41.03 

-37.53 

-33.97 

-30.36 


-50.00 166.69 -50.0593 .059255 

-40.00 172.55 -40.0605 .060527 


180 

-26.69 

-22.97 

-19.18 

-15.34 

-11.44 


- 30 .<30 170.18 -30.0336 .033602 

-20.00 183.57 -20.0041 .004107 


190 

-7.49 

-3.47 

.61 

4.75 

8.95 


-10.00 188.75 -9.9694 -.030637 

.00 193.73 .0513 -.051264 


200 

13.22 

17.55 

21.94 

26.39 

30.91 


10.00 190.52 10.0619 -.061867 

20.00 203.15 20.0596 -.059597 


210 

35.50 

40.15 

44.87 

49.66 

54.51 

— 

30.00 207.62 30.0397 -.039662 

40.00 211.94 40.0054 -.005437 

* 

220 

59.43 

64.43 

69.49 

74.62 

79.82 


50.00 216.13 49.9684 .031590 

60.00 220.20 59.9328 .067192 


230 

85.10 

90.44 

95.86 

101.36 

106.92 



70.00 224.17 69.9159 .084074 

BO. 00 228.04 79.9327 .067263 

240 

112.57 

116.28 

124.08 

129.95 

135.89 


90-00 231.83 89.9904 .009581 

100.00 235.55 100.1195 -.119519 

250 

141.92 

148.02 

154.20 

160.46 




DEG C 


E -622 


PLS SEN T 
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FRAME 16 



E -622 


PLS SEN T 
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i 

0 

2 

4 

e 

8 





0 

-79.66 

-70.25 

-76. B4 

- 75.43 

-74.03 

- : r& i al no 

MJS-2 









' ' ' r’ ONC T ] ON 

PLE modulator temp 



30 

-72 62 

-71.21 

-69 BO 

-68 39 

66.99 

F OS TREESWJ TCM ID . . 

se 



20 

-65. SB 

-64.17 

-62 76 

-61 .36 

-59.95 

FDS P*nCE 

750-700 G~MS 









fos serial no 

2 



*0 

- 58 . 54 

-57.13 

-55.72 

-54 . 32 

-52.91 

S'JE-iiv R£F NO 

20 ?c 



4 j 

-51.50 

-50.09 

-4B.69 

-47.28 

45 87 

L jz - * ■ IEP1AI NO... 

002 













so 

-44.46 

-43.05 

-4J 65 

-40.24 

-38.83 

ME *3 CAL 16 RANCE . . - 

-70 TO 100 DEO C 









MCAS CAL IB TEMP .... 

- 




-37.42 

-36.02 

-34.61 

-33.20 

-31.79 

VD K'ZZtt serial no. . . . 

LASS 



70 

-30.30 

-28.50 

-27.57 

-26.16 

-24.75 

xD- TER IMPEDANCE ... 

- 









*D. :EP CALlB DATE . . 

1 1 NOV 74 



eo 

-23.34 

-23 .94 

-20.53 

-19. 12 

-17.71 

D-*T«- FPEPAPEC Br . . . 

ROSE MONT ENG CO 


D 

A 

30 

-16.31 

-14.90 

-13.49 

- 12 . oe 

-10.67 

APPROVED Er . . . 

- 


A 







APPROVED b» ... 

P. C OTTER 


T 

100 

-9.27 

-7.86 

-6.45 

-5.04 

-3.64 

Rur. DATE 

0?97 F.C22PLS 


A 











no 

-2.23 

-.82 

.59 

2.00 

3.40 

D DEVIATION . 

. 3495022c* 00 












N 

120 

4.81 

6.22 

7.63 

9.03 

10.44 

C IE? r : C I ENT AO 

. 79657 113-02 


lj 







CGE-FiClENT 

. 70390205*00 


M 

130 

11.85 

13.26 

14.67 

16.07 

17.48 

ME-i --PAND V.'JTK-T 

Ft /) 


B 

140 

18.69 

20.30 

21.70 

23.11 

24.52 

E j On 

PPII-C 

DELTA 

F 










L 

R 

150 

25.93 

27.34 

28.74 

• 30.15 

31.56 

-*O.CC 12.74 

- 70 .6 907 

.690693 







-60.00 27. ii 

-60.4268 

.4?8842 


160 

32.97 

34.38 

35.^0 

37.19 

38.60 

-SO. GO 4J.M 

-SO. 2076 

.207632 








-40.00 56 . ;j 

-40.0271 

.027062 


170 

40.01 

41.4) 

42. B2 

44 23 

45.64 

• •**** 70 r ; 

-29.6921 

.117947 








-20.00 es Or 

-19.7675 

.232475 


180 

47.05 

40.45 

49.06 

51.27 

52.60 

- 10.00 v9>.\ 

-9 . t?96 

.311442 








.03 113.es 

. 3S99 

. 359927 


190 

54 . 08 

55.49 

56.90 

58.31 

59.72 

iO.OO W 127.92 

10.3630 

. 393015 








2 0.00 Me . ie 

20 . i? 07 

. 360699 


200 

61 12 

62.53 

63.94 

65.35 

66.75 

•o.oo ;sf.2B 

30.3*79 

. 347901 








40.00 170.40 

40.2997 

.269705 


210 

68. 16 

69 57 

70.98 

72.39 

73.79 

! 0 . 00 I94.4d 

SO. 2010 

.201029 








-3.00 196 . S3 

60. OS 70 

. 096950 


220 

75 20 

76.61 

7e.02 

79.42 

80.83 

1 . Ou 21 2 . * 3 

E9.94J4 

.057609 , 








? 3.00 226.4b 

79.7724 

. 22 r 567 


230 

62.24 

83.65 

85.06 

86.46 

87.87 

9 .‘.00 2*0.42 

09.5771 

. 422926 








ICC. -JO 2S4 31 

99 3512 

. 646767 


640 

59.28 

90.69 

92.09 

93.50 

94.91 





250 

_-6 . 32 

97.73 

99.13 

100.54 
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DEG C 
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E-623 PLS MOD T 






Jc 3 I 
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PUN 0397 0n?2PLS 


MJS77 SERIAL NO MJS-2 

«JS77 FACTION PLS ELECTRON] CS TEMP 

FDS TPEESwITCn ID . . F6 

FDS RANCE 350-700 OhTtS 

FDS SERIAL NO 2 


SUB ASS r REF NO 2032 

SU6ASSY SERIAL NO... 002 

«EAS CAL IS RAnCE ... -70 TO 100 DEC C 
MEAS CAL IB TEMP 

XDOCER SERIAL NO 6135 * 

XDUCER IMPEDANCE ... 


XDUCER CAL IB 

date .. 

3 MAR 72 


DATA PREPARED 

BY ... 

ROSE MONT ENC CO 


APPROVED 

BY . . . 

- 


APPROVED BY . . . 

R. ROTTER 


RUN DATE 
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APPENDIX C 



cpresp.f 


cpresp.f 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


C 


c 


c 

5 


20 

C 

C 


03 121 1 84 16:14:04 03/21 1 84 14:46:26 $380300 MJS 

FUNCTION CPRESP(V,UU,SUP,ZCU,AEW,III) CPRESP 

11/12/78 VMV 

REVISED 12/30 f 78 VMV 

COMPUTES RESPONSE FUNCTION OF VOYAGER CUP, IN (CURRENT UNIT)CM**3 / SC 
V(3) -VELOCITY VECTOR , CM /SEC 

UU -EFFECTIVE MODULATOR VOLTAGE (VOLTS*Z , INCLUDING SIGN) 

SUP - SUPPRESSOR VOLTAGE 
ZCU- PARTICLE CHARGE 

AEW -PARTICLE MASS (AEW*V* *2 -ENERGY IN EV) 

III-CUP INDEX 1-A , 2-B, 3-C, 4-D PROTON MODE, 5-D ELECTRON MODE 


SUBPROGRAMS CALLED TRANSNJRANSDLIPLNLISPLDARMAINARSIDE 

x - 5 


REAL* 8 $$$ DT(6) / '$$DATE$$', '03 i 21/ 84 Y 14:46:26 \ 
' CPRESP V FORTRAN',' MJS '/ 


DIMENSION V(3) 

common /comellos/ ilossm,ilossd,ilossi 
ilossm 1 for el voltage, 2 for e2 or ion voltages 
ilossd 0 no losses, 1 isotropic secondaries, 2 cosine law 
data dml / 1.6153/, ame /3.553513e-19/ 

VZ«V(3) 

IF(VZ.LE.O.)GO TO 99 
U-AEW*VZ**2 

GET TRANSPARENCY AND BEAM DISPLACEMENT FACTOR 
IF(III.GE.4)GO TO 5 
TR«TRANSN(V,U,UU,SUP) 

BDSPL-DISPLN(U,UU,SUP) 

GET AREA 

X-BDSPL*V(1)/VZ 
Y«BDSPL*V(2)/VZ 
AR«ARMAIN(X,Y) 
alossd » 0.0 
GO TO 20 

COME HERE FOR SIDE SENSOR 
TR-TRANSD(V,U f UU,SUP) 

CALL DISPII)(Vjy j UU,SUP,X,Y) 

ARsAJRSIDE^.Yi^^iilossi) 
pK< ar*aIoss^ 


.> F $ © 

:(p ©) f 



A 

Wi) 

v a') - o.o 1 

V (x) * s W 4? / 

i 

l/a)- CoS © / 

L— — 1 


FoK J'OCQ^/'OG tnotcajs 


AT JaJQlE & 


s c t I 7 u~T£7o~ 1 

I lllxI 


1 L 
(L 


osi/ - 7 

opr.,!*** 


\ i 


alossd * 0.0 


else 


\ 


\ 



CT/I/^S CO^v-:r 

A*JT> ? - h **6 *UT *' To 

Of C * l l£ CT O K. y 

/ LL UsY* O 

rt*n~ f OT 





^end If x 

(ilossd ,lfc\0) then 


vt2 - v(l)**2 + vQ)**2 
vt * sqrt(vt2) x 

vz2 » vz**2\ ^r /4 & ' 

w2 - vt2 + V^2 ' (/ 

vv - sqrt(vv2) \ a 

epem - ame*(vt2 + vz2) / t J - 

costhp - \z/ vv j c 

del\- dml rcosthp*epem**-0.35* ^ 

\ ( 1 .0-exp(-2.28 *costhp*epem* * 1 . 35)) f , , t _ { 

alossd V 1.0*del*aloss 

tL 




7i.o^s £> -- ± , 2 

l LO ^ 


ja r 


3 


/=y/£ /A/ C/2>|1>V7- 


ar * 0.0 
alossd^* -alossd 


£9 output 00/9^ r/ry /S 

(nio^Lrz) 


CONTINUE 

PUT IT ALL TOGETHER 

CPRESP-ZCU*VZ*TR*(AR - alossd) 

CHANGE SIGN OF CURRENT IN ELECTRON MODE, TO KEEP OUTPUT SIGNAL POSTT 


May 27 16:04 1987 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Page 1 of cprespf 





cpresp.f 


cpresp^f 


c 

IF(III.EQ.5)CPRESP--CPRESP 
RETURN 
99 CPRESP-O. 

RETURN 

END 


IVE 
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arside.f 


arside.f 


C 03122185 14:01:55 03/22/85 13:52:09 $380300 MJS 

FUNCTION ARSIDE(DAS,DSC, aloss, ils) 

C 12130/78 

C COMPUTES AREA JN SQ.CM, OF VOYAGER SIDE SENSOR 

C DAS -RELATIVE DISPLACEMENT OF APERTURE AND SHIELD CIRCLE CELTERS, CM 

C DSC -RELATIVE DISPLACEMENT OF SHIELD AND COLLECTOR CIRCLE CENTERS, CM 

c aloss is loss function from secondary electrons on collector 
C 

c 

ccc 

ccccccccc 

c 

C HERR SULLIVAN SOIL DIESES PROGRAM SCHREIBEN 

c REAL*8 $$$DT(6) /'$$DATE$$'/03 122/ 85' .'13:52:09', 

c | ' ARSIDE'/ FORTRAN' /MJS '/ 

COMMON / SDTABL / TABLSD(25,25),NDAS,NDSC,DDAS,DDSC 
C LOOK-UP TABLE FOR AREA OF SIDE SENSOR 
C TABLSD(IASJSC) - AREA OF SIDE SENSOR AT (DAS DSC) 

C NDASNDSC SITE OF ARRAY 
C DDASDDSC RESPECTIVE STEP SIZEA 

common /sidloss/ tablls(25,25,8) 
c label of losses for secondary electrons 

c same indicies as for tablsd 

c third index is for type of configuration 

c i-0 noloss, i-odd normal sur pressor, even reverse supressos 

c i< 5 isotropic secondary production , i>4 cosine law secondary 

c i- 1,2/5, 6 el voltages, i-3,4,7,8 e2 or ion voltages 

C 

LOGICAL INIT 
data INIT /.TRUE./ 

IF (INIT) CALL REDART 
INIT - .FALSE. 

ARSIDE - 0 

XAS - DAS / DDAS+1.0 

XSC - DSC/DDSC+1.0 

NAS - XAS 

NSC » XSC 

IF(NAS+ 1 .GE.NDAS.OR.NSC+I .GE.NDSC) RETURN 
IF(NAS.GT.O.AND.NSC.GT.O) GO TO 10 
RETURN 
10 CONTINUE 

C FOUR POINT BIVARIATE INTERPOLATION 
DXAS - XAS-NAS 
CXAS = 1.0-DXAS 
DXSC - XSC-NSC 
CXSC - 1.0-DXSC 
If ( ils .ne. 0) then 

aloss = (TABLls(NAS>NSC,ils)*CXAS+TABLls(NAS+l,NSC,ils)* 

| DXAS)*CXSC+ 

1 (TABLls(NAS t NSC+Lils)*CXAS+TABLls(NAS+l,NSC+l,iis)*DXAS)* 

| DXSC 

else 

aloss » 0.0 

end if 

ARSIDE - (TABLSD(NAS,NSC)*CXAS+TABLSD(NAS+1 > NSC)*DXAS)*CXSC+ 

1 (TABLSD(NAS,NSC+ 1 ) *CX AS+TABLSD(NAS+ 1 ,NSC+ 1 ) *DXAS)*D XSC 

RETURN 
END 


ARSIDE 
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redart.f 


redart.f 


C 01 / 16 / 85 16:04:22 01/ 16/85 10:11:39 $380300 MJS 

SUBROUTINE REDART 

C CMS VERSION OF SET ART FOR READING IN AREA TABLES PREVIOUSLY GENERATED 
C WITH ARGEN. 

C 

C FT21F001 FOR MAIN SENSOR 
C FT22F001 FOR SIDE SENSOR 
C 

C F ILED EF FT21F001 X DSN CSR $380300 ARMAIN DATA < PERM RECFM VBS 
C BLKSIZE 19069 DSORG PS 

C F 1LED EF FT22F001 X DSN CSR $380300 ARSIDE DATA (PERM RECFM VBS 
C BLKSIZE 19069 DSORG PS 
C 

c REAL*8 $$$DT(6) !'$$DATE$$' ,'01 / 16/ 85' / 10:11:39' , 

c | ' REDART',' FORTRAN',' MJS '/ 

COMMON / MNTABL/TABLMN(32,32),NSTPX,NSrPY,DSTPX,DSTPY ( XMAX,YMIN, 

* YMAX 

C LOOK-UP TABLE FOR AREA OF MAIN SENSOR 
C TABLMN(IXJY ) - AREA OF MAIN SENSOR AT X-VX/VZ 8 Y-VY/VZ 

C NSTPX, NSTPY SIZE OF ARRAY IN X,Y 

C DSTPX- STEP SIZE IN X, DSTPY- STEP SIZE IN Y 

C XMAX.YMIN.YMAX LIMITS BEYOND WHICH AREA IS ZERO 
C PRESUME XMIN-O. AREA SAME FOR +X.-X 
CCC 

COMMON / SDT ABL / TABLSD(25,25),NDAS,NDSC,DDAS,DDSC 
C LOOK-UP TABLE FOR AREA OF SIDE SENSOR 
C TABLSD(IASJSC) - AREA OF SIDE SENSOR AT (DAS, DSC) 

C NDASNDSC SIZE OF ARRAY 
C DDASDDSC RESPECTIVE STEP SIZEA 


c 


c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 


common /sidloss/ tablls(25,25,8) 
real tl (25,25),l2(25 > 25),t3(25,25) > t4(25,25) 
real t5(25 # 25), 16(25,25), 17(25,25), t8(25,25) 
equivalence (tl(l, l),tablls( 1,1,1)) 

equivalence (t2(l,l),tablls(l,l,2)) 
equivalence (t3(l,l),tablls(l,l f 3)) 
equivalence (t4(l,l),tabUs(l,l,4)) 
equivalence (l5(l,l) l tablls(l,l,5)) 
equivalence (t6(l,l),tablls(l,l,6» 
equivalence (t7(l,l),tablls(l,l,7)) 
equivalence (t8(l,l),tablls(l,l,8)) 
label of losses for secondary electrons 
same indicies as for tablsd 
third index is for type of configuration 


i»0 noloss , i-odd normal surpressor, even reverse supressos 
i< 5 isotropic secondary production, i>4 cosine law secondary 
i-l,2J,6 el voltages, i-3,4,7,8 e2 or ion voltages 


LOGICAL FIRST 
integer ia,ib 

data FIRST /.FALSE./ 
IF (FIRST) RETURN 


FIRST * .TRUE. 

open(2J file-" redart. data" form- " unformatted ') 


~ r 1 — v — 

ia - -31 



call init( 1, ia, ib) 



READ (31) TABLMN, 


DSTPX,DSTPY,XMAX,YMIN,YMAX 

1024 1 

1 

11111 

1024 1 

2 

1 2 3 4 5 

READ (31) TABLSD, 


DDAS.DD SC 
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redart.f 


redart.f 


C 

C 


c 


100 


625 111 
625 1 2 1 

read (31) 

| NDAS.NDSC 

| , NSTPX.NSTPY 

call initex( 1, ia, ib) 
read (31, end- 100) tl 
read (31, ends 100) t2 
read (31, end* 100) 13 
read (31, end- 100) t4 
read (31, end- 100) t5 
read (31, end- 100) t6 
read (31, end- 100) t7 
read (31, end- 100) t8 
return 
continue 

write (6,*) ' do loss table read y 
RETURN 
END 


1 

2 
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S/C 31 
YEAR 

DOY 

TYPE OF 
MANEUVER 

0 360 
YAW 
TURNS 

N 360 
ROLL 
TURNS 

Comment* 

— 

1978 

048 

CRSMVR 

10 

23 

ABORTED 


1978 

237 

CRSMVR 

10 

25 



1979 

034 

CRSMVR 

10 

25 



1979 

288/9 

CRSMVR 

10 

23 


- 

1988 

004 

MAG-ROLL 


1 



1988 

031 

CRSMVR 

10 

25 



1988 

121 

MAG-ROLL 

0 

1 



1988 

168 

CRSMVR 

10 

25 



1988 

189 

MAG-ROLL 

0 

1 



1988 

233 

Mini -CRSMVR 

4 

4 


— 

1980 

334 

MAG-ROLL 

0 

i 



1981 

041 

CRSMVR 

10 

25 



1981 

1 17 

MAG-ROLL 

0 

1 



1981 

134 

CRSMVR 

10 

25 



1981 

195 

MAG-ROLL 

0 

1 



1981 

321 

CRSMVR 

10 

25 



1981 

351 

MAG-ROLL 

0 

2 



1982 

089 

CRSMVR 

10 

25 


— 

1982 

218 

MAG-ROLL 

0 

2 



1982 

300 

Mini -CRSMVR 

4 

4 



1983 

1 17 

CRSMVR 

10 

25 



1983 

145 

MAG-ROLL 

0 

4 



1983 

306 

Mini -CRSMMR 

4 

4 


_ 

1984 

0 1 1 

Mini -CRSM'JR 

4 

4 



1984 

1 16 

Mini -CRSML>R 

4 

4 



1984 

193 

Mini -CRSMVR 

4 

4 


. — 

1984 

312 

Mini -CRSMVR 

4 

4 



1983 

059 

Mi ni -CRSMUR 

0 

0 

CANCELLED 


1985 

213 

MAG-ROLL tc ABCAL 

0 

4 



1983 

238 

MAG-YAW 

4 

0 



1985 

283 

Mini -CRSMUR 

4 

4 



1986 

093/4 

Mini-CRSWR 

4 

4 



1986 

153/6 

Mini -CRSM/R 

4 

4 



1986 

264/7 

Mini -CRSMUR 

4 

4 




1986 

309 

Mini -CRShf^R 

4 

4 



1987 

0 16 

Mini-CRSMVR 

4 

4 



1987 

155/6 

Mini-CRSML>R 

4 

4 

SCHEDULED 




M 366 M 366 


S/C 32 


TYPE OF 

YAW 

ROLL 


YEAR 

DOY 

MANEUVER 

TURNS 

TURNS 

Comment* 

1977 

362 

CRSMVR 

10 

25 


1978 

275 

CRSMVR 

10 

25 


1979 

155 

Mini -CRSMVR 

4 

4 


1988 

024 

CRSMVR 

10 

25 


1980 

136 

CRSWR 

10 

25 


1980 

192 

MAG-ROLL 

0 

4 


1980 

249 

MAG-ROLL 

0 

1 


1980 

341 

MAG-ROLL 

0 

1 


I960 

351 

Mini -CRSMVR 

4 

4 


1981 

077 

Mini -CRSMVR 

4 

4 


1981 

148 

CRSMWR 

10 

25 


1981 

302 

Mini -CRSMVR 

4 

4 


1982 

103 

CRSMVR 

10 

25 


1982 

293 

Mini -CRSMVR 

4 

4 


1983 

132 

Mini -CRSMVR 

4 

4 


1983 

188 

MAG-ROLL 

0 

4 


1983 

313 

Mini -CRSMVR 

4 

4 


1984 

020 

MAG-ROLL 

0 

4 


1984 

145 

Mini -CRSMVR 

4 

4 


1984 

206 

MAG-ROLL 

0 

4 


1984 

329 

Mini -CRSMVR 

4 

4 


1985 

036 

MAG-ROLL 

0 

0 

MISSED 

1985 

1 12 

MAG-ROLL 

0 

4 


1985 

150 

Mini -CRSMVR 

4 

4 


1985 

170 

MAG-ROLL t* ABCAL 0 

4 


1985 

231 

MAG-ROLL hi ABCAL 0 

4 


1985 

330 

Mini -CRSMVR 

4 

4 


1986 

036 

Mini -CRShWR 

4 

4 


1986 

1 13 

Mini -CRSMVR 

4 

4 


1986 

197 

Mini -CRSMVR 

4 

4 


1986 

273 

Mini -CRSMVR 

4 

4 


1986 

339 

Mini -CRSMVR 

4 

4 


1987 

042 

Mini -CRSMVR 

4 

4 


1987 

133 

Mini -CRSMVR 

4 

4 




MASSACHUSETTS INSTITUTE OF TECHNOLOGY 

CENTER FOR SPACE RESEARCH 

CAMBRIDOE, MASSACHUSETTS Oil }9 


June 14, 1978 

MEMORANDUM 

To: J. Belcher, G. Gordon, A. Lazarus 

From: J. D. Sullivan 

Subject: Nominal Cruise Science Maneuver (CRSMVR) 

1. A science maneuver comprises 10 yaw turns and 26 roll turns; it is 
possible to do an intial partial roll before commencing the yaw turns. We 
have requested partial rolls to bring the side sensor across the sun, and to 
sweep the main cluster through the ecliptic on different CRSMVR. 

2. The schedule for yaw turns is: 


Turn 1 

AJ1 5 

normal plasma 

2 

AJ1 5 

normal plasma 

3 

AJ03 

plasma, reversed grid 

4 

AJ1 5 

normal plasma 

5 

AJ1 0 

current cal. mod off 

6 

AJ11 

current cal. mod on 

7 

AJ1 3 

IDC mod on normal grid 

8 

AJ1 5 

normal plasma 

9 

A J 1 5 

normal plasma 

10 

AJ1 5 

normal plasma 


The purpose of the calibrations is to determine the modulator pickup with no 
plasma input and to measure the net photo current. To really clean this up, 
we probably should do a turn with IDC mod on reversed grid. 

3. The schedule for roll turns is: 

Turn 1 - 25 AJ15 normal plasma 

26 AJ10 current cal. mod off 

The prime roll turn calibration goal is intermeasurement chain calibration. 
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cprlm.f 


cprlm.f 


c 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


1 


c 

5 


10 


16 


20 

C 


03/21/84 16:14:04 03/21/84 14:46:26 $380300 MJS 

FUNCTION CPRESPjy.UUSUPTCUAE^ op foK TDC CPRESP 

) m odified to compute respons to sunligto j Q/J fachjc, Sok> 

n/i 2 r 78 — vm — 

REVISED 12/30/ 78 VMV 

COMPUTES RESPONSE FUNCTION OF VOYAGER CUP , IN (CURRENT UN1T)CM**3/SC 
V(3)-VELOC!TY VECTOR, CM/SEC 

UU •EFFECTIVE MODULATOR VOLTAGE (VOLTS*Z, INCLUDING SIGN) 

SUP • SUPPRESSOR VOLTAGE 
7£Um PARTICLE CHARGE 

AEW -PARTICLE MASS <AEW*V**2 -ENERGY IN EV) 

IU-CUP INDEX 1-A, 2-B, 3-C, 4-D PROTON MODE, 5-D ELECTRON MODE 
SUBPROGRAMS CALLED T RAN SN, TRAN SD DIPLNDISPLD ARMAINARSIDE 


REAL*8 $$$DT(6) / '$$DATE$$','03 / 21 / 84' ,'14:46:26', 

| ' CPRESP Y FORTRAN',' MJS '/ 

DIMENSION V(3) 

common /comellos/ ilossm,ilossd,ilossi 
ilossm 1 for el voltage, 2 for e2 or ion voltages 
ilossd 0 no losses , 1 isotropic secondaries, 2 cosine law 
real ara(2),alossa(2) 

data dml / 1.6153 /, ame / 3.5535 13e-19/ 
data th 1 / 40. / , th2 / 80. / , dth / 0.5 / degrad / 1 .74533e-2 / 
th - thl 
continue 
v(l) a 0.0 
U - th*degrad 
v(2) - sin(tt) 
v(3) » cos(tt) 
uu - 0.0 
sup = 0.0 
zcu - 1.0 
aew ■ 1.0 
iii » 4 
VZ-V(3) 

1F(VZLE.O.)GO TO 99 
U=AEW*VZ**2 
TRaTRANSD(V t U,UU,SUP) 

CALL DISPLD(V,U,UU,SUP f X,Y) 
ilossm*2 

do 10 ilossd- 1,2 
if (ilossd .eq. 1 ) then 
ilossi - 3 

else 

ilossi « 7 

end If 

ARaARSIDEtXXaloss, ilossi) 
ara( ilossd) a ar 
alossa(ilossd) a aloss 
continue 

write ( 6,16) th, x, y, tr, ara, alossa 
for ma t(f6,2» 1 x,2f6 . 2,2x,pe 1 3.3 >2x,2e 13 .3,2x f 2e 1 3 .3) 
if ( th .gL th2) go to 20 
th * th+dth 
go to 1 
CONTINUE 

PUT IT ALL TOGETHER 
call exit 
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cprim.f 


cprlm.f 


c CPRESP-ZCU*VZ*TR*(AR - alossd) 

c C CHANGE SIGN OF CURRENT IN ELECTRON MODE, TO KEEP OUTPUT SIGNAL POSIT 

cC JV* 

c IF(IllEQJ)CPRESP—CPRESP 

c RETURN 

c 99 CPRESP-O. 

c RETURN 

END 


/«« i2 1987 
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V 

40.00 

3.87 

1.22 

0 . 639*+00 

/ 

0.497*402 

'N 

0.497*402 

0.925*401 

0.497*401 

40. SO 

3.94 

1.24 

0 . 638*400 

0.490*402 

0.490*402 

0.915*401 

0.495*401 

41.00 

4.01 

1.26 

0.637*400 

0.483*402 

0.483*402 

0.906*401 

0.493*401 

41 . SO 

4.08 

1.28 

0.636*400 

0.475*402 

0.475*402 

0.896*401 

0.490*401 

42.00 

4.15 

1.31 

0.634*400 

0.468*402 

0.468*402 

0.886*401 

0.488*401 

42. SO 

4.22 

1.33 

0 . 633*400 

0 . 460*402 

0.460*402 

0.876*401 

0.486*401 

43.00 

4.30 

1.35 

0 . 631*4-00 

0.452*402 

0.452*402 

0.865*401 

0.484*401 

43.50 

4.37 

1.38 

0.630*400 

0.444*402 

0.444*402 

0.855*401 

0.481*401 

44.00 

4.45 

1 . 40 

0.628*400 

0.436*402 

0.436*402 

0.844*401 

0.479*401 

44 . SO 

4.53 

1.43 

0.627*400 

0.428*402 

0.428*402 

0.833*401 

0.476*401 

45.00 

4.61 

1.45 

0.625*400 

0.419*402 

0.419*402 

0.820*401 

0.471*401 

4S.S0 

4.69 

1 .48 

0.623*400 

0.410*402 

0.410*402 

0.805*401 

0.465*401 

46.00 

4.77 

1 . 50 

0. 622*400 

0.400*402 

0.400*402 

0.789*401 

0.459*401 

46.50 

4.86 

1 . 53 

0.620*400 

0.389*402 

0.389*402 

0.773*401 

0.452*401 

47.00 

4.94 

1 . 56 

0.618*400 

0.379*402 

0.379*402 

0.756*401 

0.445*401 

47.50 

5.03 

1.58 

0 . 616*400 

0.369*402 

0.369*402 

0.740*401 

0.438*401 

48.00 

5.12 

1.61 

0.615*400 

0.358*402 

0.358*402 

0.723*401 

0.431*401 

48 .50 

5.21 

1.64 

0.613*400 

0.348*402 

0.348*402 

0.705*401 

0.423*401 

49.00 

5.30 

1.67 

0.611*400 

0.337*402 

0.337*402 

0.687*401 

0.416*401 

49.50 

5.39 

1.70 

0.609*400 

0.326*402 

0.326*402 

0.669*401 

0.404*401 

50.00 

5.49 

1.73 

0.607*400 

0.314*402 

0.314*402 

0.650*401 

0.399*401 

50.50 

5.59 

1.76 

0.604*400 

0.303*402 

0,303*402 

0.631*401 

0.390*401 

51.00 

5.69 

1.79 

0.602*400 

0.291*402 

0.291*402 

0.612*401 

0.381*401 

51.50 

5.79 

1.82 

0.600*400 

0.280*402 

0.280*402 

0.591*401 

0.372*401 

52.00 

5.90 

1.86 

0.598*400 

0.268*402 

0 . 268*402 

0.570*401 

0.362*401 

52.50 

6.00 

1.89 

0.595*400 

0.256*402 

0.256*402 

0 . 549*401 

0.351*401 

53.00 

6.11 

1.92 

0.593*400 

0.243*402 

0.243*402 

0.527*401 

0.340*401 

53.50 

6.23 

1.96 

0.590*400 

0.230*402 

0.230*402 

0 . 503*401 

0.327*401 

54.00 

6.34 

2.00 

0.588*400 

0.217*402 

0.217*402 

0.478*401 

0.314*401 

54.50 

6.46 

2.03 

0.585*400 

0.204*402 

0.204*402 

0.454*401 

0.301*401 

55.00 

6.58 

2.07 

0.582*400 

0.191*402 

0.191*402 

0.428*401 

0.287*401 

55.50 

6.70 

2.11 

0.580*400 

0.177*402 

0.177*402 

0.402*401 

0.272*401 

56.00 

6.83 

2.15 

0.577*400 

0.164*402 

0.164*402 

0.376*401 

0.257*401 

56.50 

6.96 

2.19 

0.574*400 

0.151*402 

0.151*402 

0.350*401 

0.242*401 

57.00 

7.10 

2.23 

0.571*400 

0.137*402 

0.137*402 

0.322*401 

0.226*401 

57.50 

7.23 

2.28 

0.567*400 

0.124*402 

0.124*402 

0.295*401 

0.209*401 

58.00 

7.37 

2.32 

0.564*400 

0.111*402 

0.111*402 

0.267*401 

0.192*401 

58.50 

7.52 

2.37 

0.561*400 

0 . 973*401 

0.973*401 

0.238*401 

0.174*401 

59.00 

7.67 

2.41 

0.557*400 

0.842*401 

0.842*401 

0.209*401 

0.155*401 

59.50 

7.82 

2.46 

0.554*400 

0.718*401 

0.718*401 

0 .181*401 

0.136*401 

60.00 

7.98 

2.51 

0.550*400 

0.594*401 

0.594*401 

0.152*401 

0.116*401 

60.50 

8.14 

2.56 

0.546*400 

0.475*401 

0.475*401 

0 . 124*401 

0.964*400 

61.00 

8.31 

2.62 

0.542*400 

0.365*401 

0.365*401 

0.964*400 

0.762*400 

61.50 

8.49 

2.67 

0.538*400 

0.255*401 

0.255*401 

0 . 689*400 

0.556*400 

62.00 

8.67 

2.73 

0.534*400 

0 . 164*401 

0.164*401 

0.450*400 

0.368*400 

62.50 

8.85 

2.79 

0.530*400 

0.828*400 

0.821*400 

0.230*400 

0.191*400 

63.00 

9.04 

2.85 

0.525*400 

0.215*400 

0.215*400 

0 . 614*-01 

0 , 525e-01 

63.50 

9.24 

2.91 

0.521*400 

0 . 679«-01 

0 .679*-01 

0 . 197*-01 

0 . 170e-01 

64.00 

9.45 

2.97 

0.516*400 

0.000*400 

0.000*400 

0.000*400 

0 . 000e+-00 

64.50 

9.66 

3.04 

0.511*400 

0.000*400 

0.000*400 

0.000*400 

0 . OOOe + QO 

65.00 

9.88 

3.11 

0.506*400 

0-000*400 

0 . 000*400 

0 . 000*400 

0 . 00Qe*QG 

65.50 

10.11 

3.18 

0.501*400 

0 . 000*400 

0.000*400 

0 . 000*400 

0 . 0Q0e+00 

66.00 

10.35 

3.26 

0.495*+00 

0 . 000*400 

0.000*400 

0 . 000*400 

0.000e400 

66.50 

10.60 

3.34 

0.490*400 

0.000*400 

0.000*400 

0 . 000*400 

0 . 000e + 00 

67.00 

10.85 

3.42 

0.484*400 

0 . 000*400 

0.000*400 

0.000*400 

0 . 000*400 

67.50 

11.12 

3.50 

0.478*400 

0 . 000*400 

0.000*400 

0.000*400 

0.000*400 

68.00 

11.40 

3.59 

0 . 471*400 

0.000*400 

0.000*400 

0 . 000*400 

0.000*400 

68.50 

11.70 

3.68 

0 . 465*400 

0 . 000*400 

0.000*400 

0.000*400 

0.000*400 

o 

o 

VO 

12.00 

3.78 

0 . 458*400 

0.000*400 

0.000*400 

0.000*400 

0.000*400 

69.50 

12.32 

3.88 

0.451*400 

0 . 000*400 

0 . 000*400 

0.000*400 

0.000*400 

70 . 00 

12.66 

3.98 

0 .444*400 

0 . 000*400 

0.000*400 

0.000*400 

0.000*400 

70 . 50 

13.01 

4 . 10 

0 . 436*400 

0 . 000*400 

0.000*400 

0 . 000*400 

0.000*400 

71.00 

13.38 

4.21 

0 . 428*400 

0 . 000*400 

0.000*400 

0 . 000*400 

0.000*400 

71 . 50 

13.77 

4.33 

0.420*400 

0.000*400 

0.000*400 

0.000*400 

0.000*400 

72.00 

14.18 

4.46 

0 .412*400 

0 . 000*400 

0.000*400 

0.000*400 

0.000*400 

72.50 

14.61 

4.60 

0.403*400 

0.000*400 

0.000*400 

0.000*400 

0.000*400 

73.00 

15.07 

4.74 

0 .393*400 

0.000*400 

0.000*400 

0.000*400 

0.000*400 

73.50 

15.55 

4 . 90 

0 .384*400 

0 . 000*400 

0.000*400 

0 . 000*400 

0 . 000*400 

74.00 

16.07 

5 . 06 

0 .374*400 

0.000*400 

0.000*400 

0.000*400 

0.000*400 

74 . 50 

16.61 

5.23 

0.363*400 

0.000*400 

0 . 000*400 

0.000*400 

0.000*400 

75 . 00 

17.20 

5.41 

0 . 352*400 

0.000*400 

0.000*400 

0.000*400 

0.000*400 

75.50 

17.82 

5.61 

0.341*400 

0.000*400 

0.000*400 

0.000*400 

0.000*400 

76.00 

18.48 

5.82 

0.329*400 

0.000*400 

0 . 000*400 

0.000*400 

0.000*400 

76.50 

19.19 

6.04 

0.317*400 

0.000*400 

0 . 000*400 

0.000*400 

0.000*400 

77.00 

19.96 

6.28 

0.304*400 

0.000*400 

0 . 000*400 

0.000*400 

0.000*400 

77.50 

20.78 

6.54 

0.291*400 

0.000*400 

0 . 000*400 

0.000*400 

0.000*400 

78.00 

21.68 

6.82 

0.277*400 

0.000*400 

0.000*400 

0.000*400 

0.000*400 

78.50 

22.65 

7.13 

0.263*400 

0.000*400 

0.000*400 

0.000*400 

0.000*400 

79.00 

23.70 

7.46 

0.248*400 

0.000*400 

0 . 000*400 

0.000*400 

0.000*400 

79.50 

24.86 

7.83 

0 .232*400 

0.000*400 

0 . 000*400 

0.000*400 

0.000*400 

80.00 

26.13 

8.23 

0.216*400 

0.000*400 

0 . 000*400 

0 . 000*400 

0 . 000*400 

80.50 

27.53 

8.67 

0 .199*400 

0.000*400 

0 . 000*400 

0 . 000*400 

0.000*400 
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APPENDIX F 



FRINGE FIELD CALCULATIONS NEAR 
NEGATIVE MODULATOR GRID SUPPORT 


From: G. S. Gordon Jr. 

Date: January 6, 1988 

The directory /usr4/rlm/firinge contains the programs for computing the hinging field in the region near 
the negative modulator grid for the side sensor; some more generally use full programs are also included. 
All of the programs are compiled and linked using make and are run by direct invocation. The user must 
set up any environmental variables by hand, although some sample files (e.g., setft[123]) exist which can be 
run with the source command. 

The most important programs are: 

elec 

ELEC calls initvolt to read in the potential field from unit ftlO. From specified initial conditions 
(stdin), ELEC then integrates an electron trajectory through the electric field to the plane containing the 
shield ring and supressor grid (the assumed potential there is zero as in the reversed supressor 
configuration). If ft 11 is defined, the trajectory is written out on ftll in a format suitable for mongo. 
Note that the voltage is normally negative and the the external units are in inches while the internal 
dimensions are in meters. Note that TRAJ is a logical part of ELEC, and that the PARAMETERS must 
match. VOLTAGE computes the E field from the read in potential field. It should be noted the the 
field at the very edge is not correct as the first equipotential contour does not lie exacdy on the metal 
conductor, i.e., at the exact boundary. 

pltxy 

PLTXY plots out the potential contours, E field lines, or the inverse of the conformal 
map(real,imaginary), depending upon the type of plot specified: 1, 2, 3, 4. PLTXY does its plotting 
with mongo coutours, and it will plot to several devices, s=sun, t=tek, l=landscape laser, p-portrait 
laser, with different contour integrval which are documented (some what) in the program. 


wfromz 

WFROMZ inverts the function ZFUNC using a variation on Newton Rapson method. ZFUNC is z (w) 
which is the solution of 

dz _ k Vw — 1 
dw 7 i wVw — c 


with k a hc~ 1/2 and c> 1, k > 0. The solution is given by 



cosh 1 


2w — c — 1 
c — 1 


— L cosh -1 
Vc 


[c + 11 w — 2c 
[c — 1] w 


( 2 ) 


In this problem k = 0.452 inches, h = 0.312 inches, and c = 2.0988. This solution is given on p.161 of 
H. Korber [Dictionary of Conformal Representations, Dover Publications, Inc., 1957] and in problem 
10.21 on p.1315 of Morse and Feshbach [Methods of Theoretical Physics, Me Graw-Hill Book Co., 
Inc., New York, 1953]. The inversion method tries using both the numerical and analytic derivatives to 
find the solution and accepts the closer iteration; however, if both iterations are farther from the solu- 
tion than the initial guess, the method checks in the immediate vicinity of the initial guess for a better 
solution. WFROMZ also changes the size of the interval both for computing the derivative and for 
what it considers the immediate vicinity. It returns when either the max iterations have been run or the 



required accuracy achieved (defined as moving less than derr on the last iteration or if the number of 
iterations is more than jmax/2 being within a specified error of z); both W and Z of the final point are 
returned. 


ACOSH compute the inverse sin, cos, sinh, cosh for complex* 16 arguments. It shold be noted that the 
cutlines used may not be what are expected as they were made for this probleem. 

cone.f 

CONE computes and stores the potential field as is used by ELEC. 


con[fl 2 ] ? 

These programs are closely related to one another. They all compute the conformal maps much as 
CONE does, only where cone and conf compute the inverse map, con and coni compute the forward 
map; coni and con2 plot out the contours. 
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